Designing dark energy afterglow experiments 
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Chameleon fields, which are scalar field dark energy candidates, can evade fifth force constraints 
by becoming massive in high-density regions. However, this property allows chameleon particles to 
be trapped inside a vacuum chamber with dense walls. Afterglow experiments constrain photon- 
coupled chameleon fields by attempting to produce and trap chameleon particles inside such a 
vacuum chamber, from which they will emit an afterglow as they regenerate photons. Here we discuss 
several theoretical and systematic effects underlying the design and analysis of the GammeV and 
CHASE afterglow experiments. We consider chameleon particle interactions with photons, Fermions, 
and other chameleon particles, as well as with macroscopic magnetic fields and matter. The afterglow 
signal in each experiment is predicted, and its sensitivity to various properties of the experimental 
apparatus is studied. Finally, we use CHASE data to exclude a wide range of photon-coupled 
chameleon dark energy models. 



I. INTRODUCTION 

Though the existence of the cosmic acceleration has 
been confirmed repeatedly, its cause remains a mystery. 
The simplest explanation, a cosmological constant A, is 
completely consistent with the data [l|-ljl, but leads to 
more questions. Why is A some 120 orders of magni- 
tude below the Planck density? If some new physics can- 
cels this energy density, then why isn't the cancellation 
complete? Several answers have been proposed (5l-[l2|. 
with the simplest among them reducing at low energies 
to a single effective scalar field tunneling among the large 
number of local minima of its potential. Beyond these 
specific models, it is worthwhile to ask whether generic 
low-energy effective theories, possibly responsible for the 
cosmic acceleration, would predict any effects detectable 
in laboratory experiments. 

A single scalar field "dark energy" is the simplest dy- 
namical generalization of A, but most "natural" models 
mediate unscreened fifth forces which have been excluded 
over a large range of scales [l^. Thus, unless these mod- 
els are prevented by symmetry [13, from coupling 
to matter, they must possess a mechanism for screening 
fifth forces locally. Chameleon theories are scalar-tensor 
theories with potentials chosen to make their effective 
masses larger in higher-density regions of the universe, 
allowing them to "hide" from fifth force constraints [l6l - 
[22j . Symmetron theories j23l - l25j screen their fifth forces 
through a restoration of symmetry at high densities, 
while Galileons [2^, [l^] have non-canonical kinetic terms 
which reduce their effective matter couplings. 

Chameleons were the first of these screened theories to 
be discovered, and are likely the best-studied. If they 
couple to photons as expected [l^, then the very ef- 
fect which enables them to evade fifth force constraints 
also allows chameleons produced through photon oscil- 
lation to be trapped inside a vacuum chamber. Photon 
regeneration from such chameleons could produce a de- 
tectable afterglow (29l - l3l| which has been constrained by 
the GammeV and GammeV-CHASE (hereafter CHASE) 



experiments [31|-|35| . The goal of the current work is to 
study the behavior of chameleon particles in afterglow 
experiments. Specific examples are based upon CHASE, 
and constraints use CHASE data, but we aim to provide 
a general discussion of the design and analysis of after- 
glow experiments applicable to future experiments of this 
form. 

Afterglow experiments rely on two effects: oscillation 
and reflection. The rate at which a chameleon particle 
passing through a classical, macroscopic magnetic field 
oscillates into a photon, and vice versa, has been com- 
puted semiclassically [sj, ■ Although the smooth vari- 
ation of the magnetic field inside an afterglow experiment 
could lead to the adiabatic suppression of oscillation, we 
show that the quantum measurement of particle content 
made by glass windows inside the magnetic region al- 
most completely mitigate this suppression. Thus the 
chameleon production rate and the photon regeneration 
rate may be computed. 

Meanwhile, the reflection of a chameleon particle from 
a homogeneous region of high density is simply a mat- 
ter of energy conservation; a particle with a given energy 
cannot enter a region of space where its effective mass 
exceeds its total energy. We consider a real solid as a 
lattice of atomic nuclei surrounded by a nearly homo- 
geneous electron cloud and show that such a solid may 
be approximated as homogeneous for the purpose of de- 
termining whether chameleons reflect. Such reflection 
allows chameleons to be "bottled" in a vacuum cham- 
ber with dense walls, where they remain until regener- 
ating photons. Moreover, we show that averaging over 
photon polarizations washes out the dependence of the 
afterglow signal on the potential-dependent chameleon- 
photon phase, which was calculated by |37| . Thus the 
predicted afterglow signal is relatively robust with re- 
spect to the chameleon potential at high densities. 

Finally, we predict the afterglow signal for CHASE us- 
ing a Monte Carlo simulation which we cross-check 
against an analytic approximation improving upon (3^ . 
This prediction is shown to be unaffected by chameleon 
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scattering from atoms inside the laboratory vacuum, ro- 
bust with respect to surface roughness in the chamber 
walls, and relatively insensitive to chamber properties 
such as the reflectivity of the walls and the geometry of 
the apparatus. CHASE data are analyzed using the pro- 
file likelihood method [11] . The mo del- independent con- 
straints of [3^ arc reproduced and elaborated upon, then 
extended to a wider variety of chameleon models such as 
dark energy models. We then place CHASE constraints 
in context by comparing them to other chameleon con- 
straints as well as to forecasts. 

The paper proceeds as follows. After introducing 
photon-coupled chameleon theories in Sec. [Hi we study 
chameleon particle interactions with Fermions, photons, 
and other chameleon particles in Sec. lHIl Section lTvl looks 
at chameleon particle rcfiection from barriers of mat- 
ter. The computation of [s^ is corrected and improved 
upon in Sec.|Vl which uses simple approximations for the 
chameleon initial conditions and the magnetic field. A 
more accurate model of the magnetic field, and the effects 
of windows inside the magnetic field region, arc studied 
in Sec. ED Sec. IVHI compares GammeV and CHASE. A 
Monte Carlo simulation of CHASE is used in Sec. IVHIl 
to compute decay and afterglow rates as well as to study 
the sensitivity of the afterglow signal to properties of the 
vacuum chamber. Sec. llXI presents the data analysis, dis- 
cusses systematic uncertainties, and uses CHASE data to 
constrain several chameleon field models. We conclude in 
Section IS 



II. BASICS OF CHAMELEON PHYSICS 



A. Action and effective potential 



We study photon-coupled scalar chameleon theories 
with actions S = Sc/, + S-y + S'm of the following form: 
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Here, £m is the Lagrangian density for matter fields t/i* 
moving along geodesies of the metric exp(2/3m(/)/Mpi)(7^,y. 
This is equivalent to a coupling between and the trace 
of the matter stress tensor Tfj,^. Since cf) will vary by much 
less than Mpi//3m and Mp\//3^ in cases of interest, the 
precise functional forms of the couplings to matter and 
photons are not important. Expanding exp(/3^0/Mpi) 
and exp(2/3ni0/Mpi) to linear order in the expressions 
for S-y and Sm respectively, we find the usual photon and 
matter actions plus linear couplings to the chameleon 
field. 

In the presence of an electric field E{x) and a mag- 
netic field B{x) as well as a nonrelativistic matter den- 



sity p{x) = —Tj^, these couplings give the scalar field an 
effective potential 
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excluding terms suppressed by higher powers of 13(f) /Mp\. 
The scalar field equation of motion is then 
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B. Chameleon and thin-shell effects 

Consider a static matter density p{x, t) = p{x) ^ 
\E\. The scalar equation of motion ([5|) reduces to 



V, 
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-p{x). 



(6) 



If V.cj, and its derivatives are negligible, then ([6]) is sim- 
ilar in form to the Poisson equation V^^* — p{x) / {2Mp^) 
for the gravitational potential . The requirement that 
(j> and 4* remain finite as |x| — > 00 implies that cj) — 
2f3jjiMp\'^+ constant. As the size and density of the mat- 
ter distribution sourcing (j> and 5" increases, these fields 
also grow. This regime of negligible V, </> is known as 
the linear regime of the chameleon since the equation of 
motion is linear. 

Now suppose that, as (f) grows beyond a certain point, 
< begins to increase in magnitude rapidly and 
nonlinearly with (j). Then Vct> will partially cancel the 
matter source on the right hand side of (O, slowing 
the growth of the field (f>. For a sufficiently large and 
dense object, (j> will approach its bulk solution defined by 
^^('^buik) + l^mp/Mpi = 0, turning off the source alto- 
gether. This is known as the nonlinear regime; nonlin- 
earities in the equation of motion ([6]) are essential to de- 
termining the behavior of the field. The nonlinear regime 
is characterized by a rapid growth in the effective mass 

= y^l"^ of the field. Since a large mass decreases 
the range of the force mediated by i^, the field is able 
to "hide" from fifth force constraints, an effect known as 
the chameleon effect. 

Next, consider an object of constant density po- The 
change in the field from its background value, <\)^x) — 
0OO, will be approximately 2/3inAfpi(\E'(a?) — ^'oo) in the 
linear regime and will saturate at (/)buik(po) ~ 0oo in the 
nonlinear regime. We know that the nonlinear regime 
has been reached at a point x inside the object when 



|0bulk(po) - 0oo| < 2/3n,A/pi|1'(f) - 



(7) 



Since the source on the right hand side of © vanishes 
in the nonlinear regime, the field outside the object is 
effectively sourced only by the portion of the object that 
is in the linear regime. That is, the field "sees" only a 
thin outer shell of a sufficiently large and dense object. 
Such an object is said to be in the nonlinear regime or to 
have a thin shell. 
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C. Power law and dark energy potentials 

In order to provide concrete examples, we study 
chameleon potentials of the power law form 



(8) 



with (7 > and cither n < or n > 2. Given our 
convention /3„i > 0, the sign cr„ — signum(0) is pos- 
itive for n < and negative for n > 2. Eq. ([5]) is 
useful as a large-field approximation to the potential 
M|exp(KM^/(/)^) fa Ml + kM]i+^ cj)-'^ frequently used 
in the literature, in which Ma = 2.4 x 10^'^ cV is the 
dark energy scale and k is a dimensionless constant. We 
refer to the power-law-plus-constant model 



( A/4! for n = 4, 
\ kM^-" for n^i, 



(9) 



where k and A are dimensionless constants, as 
"chameleon dark energy." We will assume k = 1 for 
n 7^ 4 unless stated otherwise; such a model uses the 
same energy scale M\ in the constant and 0-dependent 
terms of the potential. 

The field which minimizes the effective potential @ in 
a bulk of matter density pmat and electromagnetic field 
Lagrangian density pc.EM = kl^-yMp^^ {\B\'^ — \E\'^) is 



'bulk = cr„ 
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Since Pmat ^ Pc,em hi almost all cases, (/)buik is typi- 
cally a function of pmat alone. Increasing pmat causes the 
magnitude of (j) to decrease for n < and to increase for 
n > 2. 

Differentiating V, we obtain ~ ncr„(7|0|"~^, which 
is always negative as appropriate to a chameleon theory. 



The effective mass Tri^a 
entiating once again: 



bulk) is found by diffcr- 



moff(0buik) = \/n{n - l)g 
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This increases with density ii n < Q or n > 2. 



(11) 



D. Chameleon-photon oscillation 

Consider the passage of a chameleon particle through 
a region of constant matter density pmat and external 
magnetic field Bq, with no electric field. Variation of the 
action ^ with respect to the electromagnetic field leads 
to 



9^ exp 
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(12) 



with the other two of Maxwell's equations unchanged. 
Oscillation between chameleon particles and photons can 



be described by perturbing about the background fields. 
Writing these perturbations in terms of a dimension- 
less chameleon amplitude ipij, and a dimensionless pho- 
ton amplitude ■0^ in the same direction as the mag- 
netic field perturbation, we have {—d^/dt^ — k'^)tpj = 
kPyBoMp^k X [ft X A;)?/'^ for a plane wave with momen- 
tum k = kk. We have assumed without loss of gener- 
ality that Bq = BqX. Here, k and x are unit vectors. 
In the relativistic, weak-mixing approximation, this is 
solved by [13, [1^ 
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Here, V- 
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(15) 
(16) 

i(t) is the probability that a particle begin- 
ning in a pure chameleon state '0-y(O) = at time t = Q 
will be a photon when measured at time t. C is some- 
times written in terms of the "mixing angle" w, with 
C = sin(2n7). Note that |ap is zero for a plane wave 
travelling parallel to the magnetic field and one for a 
wave perpendicular to the field. In the low-mass limit, 
for a wave with /c _L Sq, the oscillation probability sim- 
plifies to V-y^cf ~ PyB^f^/iiM^i). The probability for a 
particle beginning in a pure photon state to oscillate into 
a chameleon is also given by p4)) . 



E. An idealized afterglow experiment 

A chameleon particle with energy w will be excluded 
by energy conservation from a region in which its effec- 
tive mass mcff 3> (We will study this exclusion in 
greater detail in Sections IIIII and HVl ) Consider an evac- 
uated chamber whose internal mass density is pvac and 
whose walls have a density pmat- If nT-esi'PhuikiPvsLc)) ^ 
u! <^ 'Tioff(0buik(Pmat)), then the particle will propate 
freely inside the chamber, but will not be able to pen- 
etrate its walls. That is, the chameleon particle will be 
trapped inside the chamber. We refer to the condition 
u <C mefi((/!)buik(pmat)) as "chameleou containment." 

Afterglow experiments rely on oscillation and contain- 
ment to produce, trap, and detect chameleon scalar fields. 
Figure [1] shows a simple, idealized afterglow experiment. 
An evacuated cylindrical chamber has metal walls and 
two glass windows. Matter densities in the walls and 
windows are high enough that the chameleon contain- 
ment condition is satisfied. A large magnetic field Bq 
inside the chamber points in a direction perpendicular to 
the cylinder axis. 

In the production phase of the experiment, shown in 
Fig. [1] (a), photons are streamed through the chamber 
via the windows. The background magnetic field al- 
lows some of them to oscillate into chameleons. Since 
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the afterglow rate to be about four orders of magnitude 
less than the decay rate. 



III. INTERACTIONS OF CHAMELEON 
PARTICLES 



FIG. 1: An idealized afterglow experiment, (a): Produc- 
tion phase. Photons stream through chamber via entrance 
and exit windows, occaisionally oscillating into chameleon 
particles which are trapped inside. (6): Afterglow phase. 
The photon source is turned off and a detector is uncov- 
ered. Chameleon particles oscillate back into photons, which 
emerge from the chamber and reach the detector. 



these chameleon particles are trapped, a population of 
chameleons builds up inside the chamber. In the after- 
glow phase, Fig. [1] (6), the photon source is turned off 
while the magnetic field is maintained. Chameleon parti- 
cles propagating in this magnetic field oscillate back into 
photons. These regenerated photons may escape from 
the chamber through the windows, leading to an "after- 
glow" of photons from the chamber. 

In order to predict the afterglow signal expected for a 
given chameleon model in such an experiment, we must 
compute the rate Fdoc at which the chameleon population 
decays by photon regeneration, as well as the rate Fatt 
at which each chameleon particle produces detectable af- 
terglow photons. (Wc assume that photon regeneration 
is the dominant chameleon loss mode; we will consider 
another possibility in Sec. IlllPn Given these rates, the 
expected afterglow signal is shown in [34| to be 



dec 



prod^ 



-Fdoct 



(17) 



during the afterglow phase, t > 0. Here is the rate at 
which photons are streamed through the chamber during 
the production phase, and tprod is the duration of the 
production phase. 

Sections |V] and I Villi present accurate calculations of 
Fdec and Faft for realistic afterglow experiments. Here 
we can compute very rough estimates for our idealized 
experiment from Fig. [TJ Note that these estimates will 
not even be correct at the order-of-magnitude level; ex- 
perimental constraints must be based on the accurate cal- 
culations of Sec. IVIIII If the total chamber length is £tot, 
then the total chameleon-photon conversion probability 
at low TOoff , averaging over all angles, will be ^ V-yt^^ = 
/3^i3o^tQt/(4Mp[) to within a few orders of magnitude. 
The time taken for a relativistic chameleon particle to 
travel between the windows, again averaging over angles, 
will be of order £tot- The decay rate will be the conver- 
sion probability per unit time, ~ /3^i?o4ot/(4Mpj). The 
fraction of these photons reaching a detector of size rdct 



outside the chamber will be of order (rdet/^t 



The goal of this section is to compute the cross section 
for a chameleon particle to interact with other particles 
which it would encounter in the diffuse gas inside the 
vacuum chamber. In particular, we are interested in the 
chameleon-atom scattering cross section. The chameleon 
particle can interact directly with the proton and electron 
which make up a hydrogen atom. It can also scatter from 
the static chameleon field sourced by the mass density of 
the atom. We shall see that the latter effect is dominant, 
and that the cross section is approximately the square of 
the proton radius, so that chameleon-atom scattering is a 
negligible effect in an afterglow experiment. We conclude 
with a discussion of interactions between two chameleon 
particles. 



A. Scattering from Fermionic point-particles 

Consider a scalar particle of fixed mass scattering 
from a Fermion of mass mp. The matter coupling from 
([3]), with — T^' = m-pTpip, implies a Yukawa interaction be- 
tween the scalar and the Fermion with coupling constant 
gp = l^m'mF/Mpi. Assuming that /Sm ^ 10^^, this will 
be small for nucleons and lighter Fermions. Since the s- 
and t-channel Feynman diagrams have two vertices, the 
cross section <7^^^^^ will be suppressed by four powers 
of 5f- 

Chameleon-photon oscillation is strongly suppressed 
for nonrclativistic chameleons, since the chameleon and 
photon do not remain in phase. Thus we are interested 
in the case <^ p^, where is the chameleon momen- 
tum. Furthermore, chameleons in afterglow experiments 
are produced by lasers and detected in photomultiplier 
tubes, so we expect ^ 1 eV, much less than the masses 
of electrons and nucleons. The limit applicable to after- 
glow experiments is ^ pcM ^ n^F , where pcM ~ 
is the particle momentum in the centcr-of-mass frame. 
To lowest order the cross section is 



4 



247rm| ' 



(18) 



For 



1 m and rdct ^ 1 cm, we would therefore expect 



For a proton or neutron, treated as a point particle, 
this is 65I X 10^3** m^ with = 4 x lO'^^/Jm < 1- At 
fixed /?„!, (T0^_>0^ scales as m|, so the cross section for 
an electron is smaller by six orders of magnitude. Since 
electrons and nucleons are present in roughly equal num- 
bers, we may neglect the scattering between chameleons 
and electrons treated as point particles. 
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Chameleon scattering from background 
chameleon field 



Given a static matter density po{x) (assuming negligi- 
ble \B\ and \E\) with corresponding static solution (j)o{x), 
we can linearize the equation of motion ([5]) about 0o- 
With (f>{x,t) — (j>Q{x) + 5(j){x,t) wc find 



(□ - meff 5(f> = 0. 



(19) 



This tells us that, as the effective mass mcft{x) = 
™cff(0o(2?)) = Vfj,fj,{(f>o{x)y^^ varies with position, the 
total energy E = (p^ + m^^)^/^ of a chameleon particle 
remains constant. 

As a chameleon particle approaches an object with a 
thin shell, nonlinearity in V^^ causes TOeff to rise sharply. 
Energy conservation prevents a chameleon particle of en- 
ergy E from entering a region in which rries > E; such a 
particle will bounce off of the object. Thus an object of 
density po with a thin shell will scatter chameleon par- 
ticles with energies less than ?Ticff (0buik(po))- Although 
tunneling is possible, it is negligible for macroscopic ob- 
jects such as the glass windows used in afterglow experi- 
ments. 

The scattering of chameleon particles from the back- 
ground field (t>o of a massive object is used by after- 
glow experiments to trap chameleon particles. An evac- 
uated chamber with a "vacuum" of density pvac allows 
chameleon particles of energy E > mcs{4>buikiPvac))- 
These same particles will bounce from the chamber walls, 
of density pmat , if the containment condition 



E < mofr(0bulk(Pmat)) 



is satisfied. In GammeV-CHASE pvac ~ 10"^^ g/cm^ 
and pniat ^ 1 g/cm^. This difference of fourteen orders of 
magnitude means that chameleon particles with a large 
range of potentials can be trapped inside the vacuum 
chamber. 



C. Scattering from atoms 

We have studied chameleon scattering from pointlike 
Fermions as well as extended matter distributions. An 
atom is both. Here, we model the proton and the electron 
cloud as uniform-density spheres of radius rp = 0.83 fm 
and rBohr = 0.529 A, respectively, in a laboratory vac- 
uum of density pvac = 10~^^ g/cm"^. We compare the 
resulting cross sections to those obtained in Sec. IIII Al 

In our approximation, the electron cloud, with 
mass TOoiec = 511 keV/c^, has a density Poioc = 
3TOcioc/(47rr|^ijj.) 1.5 x lO^^g/^^j^s^ Defining the 
gravitational potential to be zero at infinity, we have 
«'cloc(0) = -peiecrLhr/(4M2i) = -1.9 X lO'^^ at the 
center of the cloud. The proton has mass rriprot = 
938 McV/c^, mean density pprot = 7.0 x lO^^g/cm'^, and 
gravitational potential ^'prot(O) = —2.2 x 10"'^^. The 
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FIG. 2: Minimum value of the matter coupling fim required 
for the electron cloud and the proton, respectively, to have 
thin shells, in the case of a chameleon dark energy V((/)) — 



(20) A/i +M^-"|0|". (Top): n>2. (Bottom): n < 0. 



electron cloud and the proton, respectively, will have thin 
shells only if 



|0bulk(Pcloc) - 0c 
|0bulk(Pprot) - 4>a 



I < 2/?^Mpi|*elec(0)|. (21) 
<2/3„,Mpi|*prot(0)|. (22) 



Recall that an object with a thin shell represents a large 
perturbation to the background chameleon field, from 
which incident chameleon particles may scatter. 

Consider power law potentials of the form ([8]). If n > 2, 
then the left hand sides of (f2T| [22]) will be dominated 
by |0buik|- For chameleon dark energy, V{4>) = M\ + 
M|""|0|", the thin sheU condit ion becomes 



i|MlA/pi 



/ A/a 



V2A/pi|*| 



(23) 



with the density p and gravitational potential 5* appro- 
priate to each object. Figure[2]( Top) shows the minimum 
values of (5^ necessary for each object to have a thin 
shell. If n < 0, then the left hand sides of (HH [21]) will 
be dominated by the background field value l^ool- As- 
suming that the vacuum chamber containing the atom 
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is much larger than the chameleon Compton wavelength 
Weff (0buik(Pvac))~^ at the vacmim density Pvac, the thin 
shell condition is just ()23p with p = pvac- Fig. [2] {Bot- 
tom) shows the minimum /3„i required for a thin shell in 
the n < case. 

Colliders already exclude /3„i > 10^^ for similar 
chameleon models (28l |. In the remaining parameter 
space, the electron cloud will only have a thin shell when 
n is tuned to be very close to zero, —10"^ < n < 0, 
and we do not consider this case further. The proton 
does have a thin shell in a substantial fraction of the 
allowed parameter space. Thus a chameleon particle in- 
cident upon a hydrogen atom will pass right through the 
electron cloud, which represents only a small perturba- 
tion to the background mass, and will interact solely with 
the proton. 

In order to find the background chameleon field (t>o{r) 
due to the proton, we solve ^ for the spherical tophat 
density distribution by which we approximate the pro- 
ton. Since the boundary conditions 0o(O) — and 
limr_>.oo = 4*00 are defined at different r, we solve (jH]) 
numerically using the shooting method: we guess a value 
of (/'o(O), solve to find the field at large r, and refine 
our guess. Since scattering will be important only when 
the proton has a thin shell, TOoff ('?!'buik(Pprot))?' ^ 1, we 
focus on this regime. We immediately run into numer- 
ical difficulties due to the exponential sensitivity of the 
large-r field to the central value. This sensitivity can 
be seen by linearizing (j)o{r) = </)buik(/Oprot) + S(j){r) in- 
side the proton, resulting in \7'^6(j) = m^ff prot^'^ '^^^^ 
"leff.prot = "^off((/'buik(Ppiot))- Thc hucarizcd equa- 
tion has the one-parameter family of solutions 6(j){r) = 
i^sinh(TOeff,prot'^)/(TOoff,prot'^), Valid as long as \6(j)\ < 
|</'buik(Pprot)|- Although this linear approximation will 
be inapplicable at the proton's surface, r ~ rp, it will be 
valid for r < rp — A?' for some Ar. In the thin shell case, 
we can choose ^^^^ <^ Ar <C rp , begin our numerical 
solution at rp — Ar using the linear approximation, and 
apply the shooting method to F rather than 0(0). 

Figure [3] shows the result of this calculation for a 
chameleon dark energy with n = — 1 and /3m = 10^^. For 
this model, (/ibuik (Ppiot ) = 2.5 x 10"^'' eV and moff,prot = 
9.9 X 10^® eV, implying a Compton wavelength of 
"^(Tff prot = 2.4 X 10~^rp. Choosing Ar w 10~^rp, we find 
using the shooting method that log(F/0buik(Pprot)) = 
-8.58 X 10^. At r = rp - Ar, log |<50/(/.buik(Pprot)| = 
—8.17 X 10^, well within the regime of validity of the 
linearized equation of motion. 

This numerical result is consistent with theoretical ex- 
pectations. Since the effective mass inside the proton 
deviates from TOeff.prot for r > rp — Ar, the proton has 
a thin shell of thickness < Ar, which includes a frac- 
tion sa 3Ar/rp of its mass. The surface gravitational 
potential of this mass shell is 4'shcii = — 8™a/^* r'^ ^ 
—4.5 X 10"'*^, or of order 1% of the total gravitational 
potential. Then 2/3m|*shoii| = 0.2 eV, of the order of 
4>oo = 4'hu\k{Pvac) = 0.07 cV. Therefore this shell alone is 
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FIG. 3: Chameleon field and effective mass for a proton, 
approximated as a uniform-density ball of mass mprot = 
938 MeV/c^ and radius rp = 0.83 fm. The potential V((/)) 
is that of a chameleon dark energy with n = — 1 and a matter 
coupling /3m = 10^^. A vertical dotted line shows the surface 
of the proton. 



enough to saturate the chameleon field inside the proton. 

Suppose that a chameleon particle approaches the pro- 
ton from large r with an energy ^ 1 eV char- 
acteristic of a laser oscillation experiment. Evidently 
from Fig. [3l the effective mass will rise to equal Ecj, at 
r w rp — Ar sa rp. Thus the incoming chameleon par- 
ticle will scatter off of a background chameleon "ball" 
of radius approximately rp, implying a cross section 
of 47rrp. In our coupling constant regime of interest, 
f^ni 10^^, for which ^ 1, this will be much larger 
than (70.0^00 Sp/'^prot ^ 5f^p from HH). There- 
fore, the cross section for chameleon-hydrogen scatter- 
ing is dominated by semiclassical scattering of chameleon 
particles by the background chameleon potential of the 
proton, which is well-approximated by a hard sphere of 
radius rp . This is true for any chameleon model for which 
the proton, but not the electron cloud, has a thin shell. 

For laser oscillation experiments such as GammeV- 
CHASE, chameleon-atom scattering has a negligible ef- 
fect on the chameleon afterglow signal. The chameleon- 
atom scattering rate for a vacuum density of pvac ~ 
10^14 g/cm"^ made up almost entirely of hydrogen atoms 
is an — 47r7'pPvac/'7iprot = 1-6 X 10~^IIz, or of order one 
scattering event per chameleon per 10^^ passes through 
a 10 meter chamber. Thus we expect atom scattering to 
correct our predicted afterglow signal at the 10^^^ level. 
Incidentally, even if the electron cloud also had a thin 
shell, and we replaced rp in the above expression by the 
Bohr radius, the rate would go up by (rBohr/'''p)^ ~ 10^". 
Thus atom scattering would still only be a percent-level 
effect. 



D. Chameleon fragmentation 

Fragmentation is the process by which one chameleon 
particle interacts with another to produce more than two 
chameleon particles. Since photomultiplier tubes are sen- 
sitive only to energies of order 1 eV, repeated fragmen- 
tation would result in a population of chameleon parti- 
cles with too little energy to produce detectable photons. 
Thus fragmentation lowers the signal expected in after- 
glow experiments. 

Let the cross section for chameleon fragmentation be 
ffrag- The fragmentation rate of chameleons at an inital 



energy of w, Tft-ag = iV^^Vfrag/V- 



-vac = M"'sfrag, dcpCuds 



on the number -/vj,"'' of chameleons at that energy as well 
as on the volume Vvac of the vacuum chamber; Sftag is 
the cross section per unit volume. 

Next we compute as a function of time for a 

typical afterglow experiment. At initial time — tprod, 
the number of chameleon particles is zero. During the 
production phase, — tprod < ^ < 0, photons of energy 
a; ~ 1 eV are streamed through the vacuum chamber 
at a rate F^, which is about 8.8 x 10^® Hz in CHASE. 
Each of these has a probability V-y^d) of oscillating into a 
chamelelon particle. At time t ~ 0, the photon source is 
turned off, and the population of detectable chameleons 
decreases from its peak of Njl^\o) due to photon regen- 
eration (with a rate Fdoc) as well as fragmentation. At 
t > 0, the afterglow phase of the experiment, a detector 
outside the chamber can look for evidence of chameleon- 
photon oscillation. 

The evolution of A^i'^'' is described by 



dt 



F^V-,^^Q{-t) - TdccA^i"^ - SfragA^r'" (24) 



where 9 is the step function. In the production phase 
t < 0, 



2f^7^^^^sinh( ^'°^<*;*^"^' ) 



Ftot cosh ( rto.(Wp,od) ^ + Tdcc sinh (^^ 



ot(t+tprod) 

2 



(25) 

where T^^^ = V\^^^ ^F^V^^^Si-^^^. In the afterglow 
phase t > 0; 



(26) 



At times much smaller than the decay time, < 
t <C rj^"^^, the chameleon population is approximately 

^'i"^(0)/(l + (sfragAi"^(0) + Tdcc)*); fragmentation will 

dominate over decay if SfragA^0"''(O) ^ Edoc- Since the 

fragmentation rate decreases with A^^"'' , decays via pho- 
ton regeneration will eventually come to dominate, lead- 
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FIG. 4: Fragmentation processes in (j)* chameleon models. 



ing to a chameleon population of N^\o)e ^^aat^^-^^ _|_ 

SfragiVi"^(0)/Fdcc). 

Thus far our discussion has used the purely phe- 
nomenological parameter crfrag. For the particular case 
of (/)^ theory, we can estimate the fragmentation cross 
secton from processes such as those in Figure |4] using di- 
mensional analysis, Cfrag = aftagA^/w. Here atrag <C 1 
is a dimensionlcss numerical factor resulting from an in- 
tegral over the four-body phase space of outgoing parti- 
cles. We assume afrag = 1 in order to be conservative; 
as we shall see in Section IIXI CHASE constraints are 
not competitive with Casimir force constraints for this 
potential, so a lengthy numerical calculation of afrag is 
unwarranted. This ordcr-of-magnitude calculation shows 
that fragmentation cannot be neglected in CHASE for 
A > 0.001. 



In principle, any potential of the forms (|8l9p can be ex- 
panded in Taylor series about the expectation value 0buik 
of the field in a given matter density. Consider chameleon 
dark energy with n ~ —1. The ^-th order term in the 
series for £ > 4, (-l)'M50-^(f+i) (0 - ^buik)' = Vi{^ - 
0buik)^, which describes a process in which two chameleon 
particles fragment into ^—2, should contribute a quantity 
ae ^ |«,|2w2£-io ^ (MA/c.)iO(c./0buik)'7</'buik to en,,,, 
by dimensional analysis. The prefactor {M\/uj)^'^ ^ 
lQ-30 CHASE. However, for typical parameter val- 
ues, 0buik can be less than w, implying that cr^ — > oo as 
£ — > oo. The same is true for other n < 0; ^buik is typ- 
ically within several orders of magnitude of the energy 
scale A/a in the potential. 

The underlying problem is that chameleon theories are 
low-energy effective field theories whose cutoff energies 
can be below the particle energy w ~ 1 eV in an afterglow 
experiment. A proper calculation of the fragmentation 
cross section would require a more fundamental theory, 
but UV completions of chameleon theories are not yet 
well- understood [1^, so such a calculation is beyond the 
scope of this paper. Henceforth we neglect chameleon 
fragmentation for all potentials other than ^(^^. 
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IV. REFLECTION FROM A BARRIER 



A barrier as a lattice of atoms 



100000 



10000 



Section IIII Bl showed that a chameleon particle of en- 
ergy E will be excluded by, and will therefore reflect 
from, a region in which the effective chameleon mass is 
greater than E. For a homogeneous object of density po 
which satisfies the nonlinearity condition ([7]), this reflec- 
tion condition becomes moff ((/)buik(po)) > E. 

However, ordinary matter is not homogeneous, and ho- 
mogeneity is not necessarily a good approximation in the 
case of a nonlinear field such as the chameleon. Mota and 
Shaw, in references [111, [2^ , considered matter as a cu- 
bic lattice of homogeneous spheres of radius rg, with a 
lattice spacing do ^ tq. Since the spheres approximate 
atomic nuclei, we expect ~ 1 fm and do ~ 1 A. The 
mean density po = M^dQ^ of such matter depends on the 
mass Mq of each sphere as well as the spacing do, but 
is independent of vq. Suppose that at some large value 
of ro one sphere considered individually is in the linear 
regime of a chameleon theory. As rg is decreased beyond 
a certain point, the sphere will acquire a thin shell, and 
the chameleon will effectively "see" only a fraction of it. 
Thus the nonlinear chameleon theory will be sourced by 
a density smaller than the mean density poi and its mass 
will be smaller than Woff (<^buik(po))- Refs. [2l|, [g^] ap- 
proximate the actual chameleon mass inside such matter 
as the lesser of Woff (0buik(po)) and 



(27) 



with q{n) = min(l, for a power law or chameleon 

dark energy potential. 

Well inside the nonlinear regime, when the fraction 
of each nucleus seen by the chameleon is much less than 
Wcicc/^prot ^ 10~^, the chameleon will effectively ignore 
the nuclei and be sourced by the approximately homoge- 
neous gas of electrons making up the solid. This gas has a 
density po,e = Po(weiec/"^prot)(-2o/A) where Zq and Aq 
are, respectively, the atomic number and mass number 
of the atoms making up the lattice. Even a highly non- 
linear chameleon field should have a mass no lower than 
TOeff (0buik(/3o,o))- Thus wc approximate the chameleon 
mass in such a lattice by 

miattico ~ max Woff (0buik(po,o)), 

min(mcrit, Woff (0buik(/5o))) • (28) 

Figure [5] ( Top) shows miatticc as a function of the mat- 
ter coupling /3ni for several power law dark energies. 
Three regimes are seen for each model. At low /3m the 
chameleon effectively sees a homogeneous solid, and the 
mass corresponds to the bulk field value expected for such 
a solid. At intermediate /3m the mass flattens out as a 
greater portion of each atomic nucleus is "hidden" be- 
hind a thin shell. At the largest /3m, the homogeneous 
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power law index n 

FIG. 5: Chameleon mass in matter, modelled as a lattice of 
spherical nuclei in a homogeneous gas of electrons. In each 
plot, the horizontal dotted line shows m — to — 2.33 eV, the 
chameleon energy in CHASE. ( Top) Approximate chameleon 
mass. (Bottom) Critical mass of refs. [2ll. |22|. the mass above 
which the inhomogeneous nature of the matter becomes ap- 
parent to the chameleon field. 



electron gas dominates the density seen by the chameleon 
field, causing the chameleon mass to rise with /3m once 
again. 

In each case, the mass required for chameleon contain- 
ment, meft = uj = 2.33 cV in the case of CHASE, lies in 
the low-coupling regime. Figure [S] (Bottom) generalizes 
this conclusion to chameleon dark energy models of ar- 
bitrary n. Thus we are justified in treating the chamber 
wall as a homogeneous solid of density po for the purpose 
of testing the containment condition. 



B. Phase change due to reflection 

Let there be a vacuum p = in the region x > 
and a constant density Pmat at a; < 0. This is a simple 
model for a planar slab thick enough that the chameleon 
reaches its bulk value. Given this density and a power 
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law potential, ([6]) is solved in the region a; > by 
<Ao(a;) 



rricsix) = TOoff 



^surf [l + V5/2|n - '2\(f>suri ^ 
(1 - l/n)(/)bulk(Pmat) 

' - 2|mcff,surf X 



1 



m-efF,surf 



gjn - 1) 



ji-i 



jigAIpi 



1) 



(29) 
(30) 

(31) 
(32) 



where a subseript "surf" denotes the value of a function 
at the surface x = 0. 

We wish to find solutions to ([5]) representing plane 
wave perturbations to 4'o{x) far away from the wall, 
-iu:t ^j^gjj a; 0. With the definition 



Do oc e 



(/)(a;,t) = <f)o{x) + e 



(x), reduces to 

1 OJ^ 



{x + £w/TOcff,surf)^ 



(33) 



with the length £ = oj ^^2n(n — l)/|n — 2| of order the 
wavelength of the plane wave. This is solved by 



a 



= i^/^ (ciJ„(i) +C2r„(i)) 

1 = w(a; + ^a;/moff,surf) 
^2^2,^ 1/4 



(34) 
(35) 
(36) 



where J a and are, respectively, Bessel functions of the 
first and second kind, of order a. At the surface of the 
wall, a; = 0, so i = ^a;^/meff,surf- We expect ~ 1 while 
w / TOeff^surf ^ 1 , SO i w at the wall. Thus the coefficient 
of the irregular Bessel function must be small, C2 ^ c\. 
Far from the wall, will therefore be dominated by the 
regular Bessel function: 

^<i>i&Y{x) « c\\fujxJa(yix) (X ci sin(wx — 7rQ;/2 + 7r/4). 

(37) 

Writing this perturbation as the difference of an incoming 
wave and an outgoing wave with a X^(0)-dependent phase 
shift ^y, in the relativistic limit uj = k, we have 



(50(a;) cx e 



-2ie-*«/2 sin(wx-Cv/2). 



(38) 



Equating the above to (|37)) gives the phase shift for 
potentials (|8l9p . 



TT a 

2 



3n- 2 



2n 



(39) 



which we note is independent of the energy oj of the in- 
coming particle. In a real experiment particles will ap- 
proach the wall with a range of incident angles. Suppose 
now that the incoming chameleon wave has a nonzero 
angle of incidence, and let the xy plane be the plane of 
refiection. Since the y momentum ky is unaffected by 
the refiection, we may factor out e^^^^ as well as the time 



dependence: (j){x,t) = (fiQ^x) + e^^'«y~^'^*5(t){x). Defining 
fcj. = (w^ — fc^)^/^, we see that to in ([55)1 is replaced by 
kx- Since (|39| is independent of w, we have shown that 
the phase shift for power law and chameleon dark en- 
ergy potentials is independent of incident angle as well 
as energy. 

More generally, if moff(x)^ = moff (^0(2;))^ = A{x + 
D)~^ for some constants A and D, as in (PT|) . then the 
chameleon reflection phase — 7r((yl -t- 1/4)^/^ — 1/2) 
will be independent of incident angle and energy. Since 
■j^(l)o{x) = -^yi^') in the vacuum, we have m^j^ = 

i^y^^) - = '^I'-t - A{x + D)-\ This 

will be satisfied if d4>o/dx oc {x + D)'^ for ( satisfy- 
ing C(C - 1) = ^, that is, if ^o{x) (x {x + D)^+^ or 
4>o{x) (X log{x + D) up to an additive constant. Substitu- 
tion into -^(jioix) = -£j;y{4>) then implies V{(f>) — gcj)"' or 
y{4') = 5exp((^/A'f), up to additive constants, for con- 
stants g, n, and Af. Thus we have extended our re- 
sult about the independence of on incident angle and 
energy to exponential potentials. Incidentally, choosing 
C = — 1, for an exponential potential 

V{(t>) = g exp{(t>/M) + const., (40) 

results in the phase shift 

?oxp(0) = TT. (41) 

However, for more general potentials, the phase will 
depend on the incident angle and energy. A chameleon 
particle with greater momentum k^ perpendicular to 
the wall will approach closer to the wall, allowing it to 
probe a different region of the potential V{(j)) than a 
lower-fca: particle. Although exact solutions do not ex- 
ist for general potentials, we can compute (j)o{x) numer- 
ically using -^Mx) = -^Vifp)- Given (j)o{x) we may 
find mes{x) and then solve numerically for S4>{x) using 
^54>{x) = {mefi{xY — k'^)S4>{x)- Comparison of the ze- 
ros of this numerical S(j){x) to those of ([55]) is used to 
find the phase shift ^y; far from the wall, it will con- 
verge to a constant. Fig. |6] (Top) applies this technique 
to the n = — 1 power law model, verifying the energy- 
independence of £,v shown in ([39|) . Furthermore, the 
phase values converge to ^^-i = 7r/3 from (|39p . shown 
as a dotted line. 

An interesting example is the "exponential-inverse" 



potential F(0) = Mj 



which is frequently used in 



the literature. At low energies, (p ^ M\ and the po- 
tential will approximate the n = — 1 chameleon dark en- 
ergy, eq. at high energies, it will differ from 
Fig. [6] {Bottom) shows the phase shift for several differ- 
ent chameleon energies. As expected, for uj < M\ = 
2.4 X 10^^ eV, the potential approaches the n = — 1 
chameleon dark energy, and the resulting chameleon 
phase approaches ^^-i = 7r/3. At larger energies ^y dif- 
fers substantially from 7r/3. 

In the limit that the wall in an actual afterglow ex- 
periment is smooth, the afterglow rate will be dominated 
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FIG. 6: 



Chameleon reflection phase vs. 

12 „„j „ _ 1 „/„„3 



wall, with Pm = 10^'' and pmat = 1 g/cm' 



distance from 
(Top) n = -1 

chameleon dark energy V{<j)) = Ma(1 + MA/(t>). {Bot- 
tom) potential V{4>) = M^exp{M a /(/))■ In both cases, the 
dotted line shows = tt/S, the expected phase in the x — > oo 
limit for n = — 1 chameleon dark energy. 



by chameleon particles bouncing with grazing incidence, 
k^/uj < radius / length - lO'^ ^ fc^ < IQ-^ eV. Thus 
the exponential-inverse model from Fig. [6] can reasonably 
be approximated as an n = —1 chameleon dark energy 
© for the purpose of computing the phase, although this 
approximation should worsen at larger /3m- 



For potentials inconsistent with this power law ap- 
proximation, it is necessary to compute numerically the 
chameleon reflection phase as a function of the in- 
cident angle, the coupling /3i„, and the parameters of 
the potential. Although simple, this calculation is time- 
consuming and must be repeated for each new potential. 
In the absence of a compelling reason for choosing a dif- 
ferent potential, we restrict ourselves henceforth to the 
potentials (|8I9I40|) . 
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FIG. 7; Norms and arguments of the complex amplitude re- 
flectivities Aretes and Arcf,p. Colored lines assume the param- 
eters ni — 1.6, hi = 3.2, and /vis = 0.82 appropriate to the 
CHASE chamber walls. Dotted lines show the reflectivities 
in the thin-skin case; the two amplitudes are equal and real. 
See the text for a discussion of the conventions used here. In 
brief, our S and P directions differ from the standard def- 
inition, and our phases are normalized so that = (.p at 
normal incidence. 



C. Reflection and absorption of photons 

The standard treatment of photon reflection from a 
conductive medium adds an imaginary component to the 
index of refraction of the medium, ni —ihi. Assume that 
the incident wave propagates through a medium whose 
index of refraction ng ~ 1 is purely real. For stainless 
steel like that in the CHASE chamber walls, ni = 1.6 
and fii = 3.2. These quantities imply a mean reflection 
probability of /ref = 0.65. The actual measured value is 
/rcf = 0.53, or 18% lower; the unpolished chamber walls 
are less reflective than polished stainless steel. We model 
this discrepancy by assuming that 18% of the chamber 
wall is obscured and perfectly absorbing. Thus we mul- 
tiply the computed reflection probabilities by a visibility 
correction factor /vis = 0.82. 

The complex index of refraction implies that the pho- 
ton reflection occurs over a skin depth Si = a;~^nj~^. 
This introduces a phase shift which depends on the po- 
larization direction as well as the incident angle. Fig. [7] 
plots the amplitude reflectivities A-Gf,p = |^ref,p|e'^^ 



and ^rcf.s = |^rof.s|e*^^ for the P and S polarizations, 
respectively, as a function of incident angle. (Recall that 
we keep track of the magnetic held component of the 
photon, cx SB, so that our P and S polarizations 
arc switched with respect to standard conventions.) The 
norm of the amplitude reflectivity is the square root of 
the reflection probability, and the complex phase is the 
phase shift of the photon due to reflection. In our con- 
vention, this phase shift is measured relative to that of 
a chameleon particle with Dirichlet boundary conditions 
and no additional phase shift; that is, our conventions for 
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FIG. 8: Particle in an afterglow experiment. The magnetic 
field region is shaded. Inside this region, where chameleon- 
photon oscillation takes place, the particle trajectory is drawn 
as a dashed path. Next to each dashed segment is shown the 
contribution of that segment to the total photon amplitude. 
From [H. 



the photon phase shift differ from standard conventions 
by TT. 

Reflection becomes perfect in the thin-skin limit hi — > 
oo. ^ref.p and Ai-of,s both approach unity in this 
case. Ref. [13] implicitly used this approximation, set- 
ting /vis = 0.53 to match the observed reflection prob- 
ability. Fig. [7] shows that, with this value of /vis and 
hi = 10^, the thin-skin limit of that reference is repro- 
duced. Henceforth we use the values ni — 1.6, hi = 3.2, 
and /vis = 0.82. 

We will see that this polarization-dependence of the 
photon phase makes Faft nearly independent of . The 
afterglow rate is dominated by particles bouncing at graz- 
ing incidence, for which ^5 and differ by w tt. Aver- 
aging over polarizations washes out the dependence on 
additional phases such as ^y. 



V. CHAMELEON-PHOTON OSCILLATION: 
ANALYTIC CALCULATION 

A. Decay and afterglow rates 

Scalar-photon oscillation was studied in [s^ and ap- 
plied to afterglow experiments in [sj]. Although [sj] did 
not include the photon polarization-dependent reflectiv- 
ities just described, the data analysis in [3^ did incor- 
porate this effect. Here we summarize the calculation 
of [U and add a photon polarization-dependence to the 
reflectivity and phase shift. 

Figure [8] shows a cylindrical oscillation chamber with 
windows at the ends and a magnetic field region (shaded) 
offset from the windows. The coordinate system used is 
identical to that of [s^. The origin lies at the center of 
the entrance window. The z axis is the cylinder axis of 
the vacuum chamber, with the unit vector z chosen to 
point in the direction from the entrance window to the 
exit, a; is a unit vector in the direction of the external 
magnetic field Bq, orthogonal to the cylinder axis, and 
y = z X X. 

Consider a particle beginning in a pure chameleon state 
at the center of the entrance window. This simplifies the 



probcm considerably by maintaining cylindrical symme- 
try. Ref. 3J] shows that, for m^g ^ Attuj/R, where uj 
is the particle energy and R the chamber radius, aver- 
aging over initial positions on the entrance window may 
be neglected for the purpose of computing the afterglow 
rate. Although a particle beginning at the center of the 
entrance window is more likely to reach the detector out- 
side the chamber, this is a purely geometrical effect for 
which we can correct quite easily. 

Since the total momentum of the chameleon particle 
is k, there is a two-parameter family of initial conditions 
beginning at the center of the entrance window. We pa- 
rameterize the initial conditions by the angle 9 between 
k and the z axis, cos{6) = k ■ z, and the angle ip between 
the X axis and the projection of k onto the xy plane, 
cos((p) = k ■ x/\k — (k ■ z)z\. (Recall that the x and z 
axes are specified by the background magnetic field and 
the cylinder axis, respectively.) Let and Lp be specified 
for a particle. Then the total number N[9) of bounces 
inside the chamber, the first bounce ni,{9) inside the re- 
gion with magnetic field i?o, and the final bounce 71r.(0) 
in the Bq region depend on 6 but not Lp by symmetry. For 
the path shown in Fig. [SI tt-l = 2, riR = 5, and iV = 6. 
As shown in ref. [s^], the photon amplitude at the exit 
of the Bq region is found by summing up the contribu- 
tions from each segment inside the Bq region shown in 
Fig. |S1 properly accounting for photon absorption and 
chameleon-photon phase shifting at each bounce. 

Eq. (fT51) gives the contribution of each segment to 
the photon amplitude. Up to a phase factor, this is 
C sin(TO^gi/(4fc))a(fc), with t the time spent inside the 
Bq region and a expressed in the basis of S and P polar- 
ization states; recall that a and C are defined by p^[T6|) . 
Since each segment in Fig. [5] beginning and ending in- 
side the Bq region has the same value of t by symmetry, 
only three values of t are necessary: the time taken by 
the particle between Bq region entry and bounce til ; the 
time tyi between two bounces inside the Bq region; the 
time tft between bounce and the Bq region exit. We 



define sl 
that 



sin(tL), sm — sin(iM), and sr, = sin(tR,), so 



mlaiZn^ -4) 



4fc cos ( 



2k sin 9 



sm 



Vff 



{ii + L- z„„ ) 
4fc cos 9 



(42) 
(43) 
(44) 



where Zn is the z value of the nth bounce, £1 is the dis- 
tance from the entrance window to the beginning of the 
Bq region, £2 is the distance from the end of the Bq re- 
gion to the exit window, and L is the length of the Bq 
region. 

Chameleon-photon phase shifts come from three dif- 
ferent sources: a potential-dependent chameleon phase 
shift due to reflection from the walls, computed in 
Sec. IIVBI a photon polarization-dependent phase shift 
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due to a nonzero skin depth during reflection, com- 
puted in Sec. IIV CI and a phase shift during propagation 
?propagation(i) through the Bq region due to a chameleon- 
photon mass difference, computed in [3J| and included in 
the time-dependent phase factor in (fT3)) . 'Cpropagation(0 
must be computed over the time intervals iL, ^m, and % 
defined earlier, leading to 



Cm 
Cr 



ff 



4A; cos 9 



UL - -] Az- li 



m 



k sin 9 



4fc cos 9 



Az - 4 - i 



(45) 

(46) 
(47) 



where Az = 2Rcoi9 is the z distance between bounces. 
Following the convention of [33 | , we redefine the complex 
reflectivities to include and the chameleon phase shift 



Ap = Aref,p exp(iCM + i^v) 
= Acf,s exp(zCM + iCv)- 



(48) 
(49) 



Next, we compute the components of a(fc) in the ba- 
sis of S and P polarization states, a and k are constant 
between bounces. Let fc„ be the particle direction before 
the nth bounce, and let a„ = a(kn)', thus a subscript 1 
denotes initial values. We have ai = (1— sin^ 9 cos^ 
{siv? 9 sin(/3 cos (p)y—(s\n9 cos 9 cos ip)z. A bounce from 
the chamber wall switches the sign of the z compo- 
nent of a while leaving the x and y components un- 
changed: (22 = fli — 2{di ■ z)z. Thus a„ will equal ai 
for odd n and 0,2 for even n. Since the plane of inci- 
dence for each bounce is spanned by z and ki , we define 
the S polarization direction by S = (z x fci)/|z x fci|; 
this remains the same during all bounces. For each n, 
the set (5, P„, k) must be orthonormal, with P„ the 
P polarization direction before the nth bounce. We 
choose the sign of P„ such that S* x P„ = fc„. With 
these definitions we find S = ~{svinp)x -\- {cos ip)y and 
Pi ~ ~{cos9 cos Lp)x — (cos 6* sin (p)y + {sm9)z; P2 has 
the same x and y components as Pi, but its z compo- 
nent has the opposite sign. Finally, since S has no z 
component, and since the z components of a and P both 
change after a bounce, we have 05 = a„ • S' = — simp 
and ap = a„ • P„ = — cos cos (p independent of n. Thus 
a„ = + apP„. 

We proceed to compute the photon amplitude at the 
exit of the Bq region step by step. Immediately be- 
fore the first bounce in the Bq region, nL, the pho- 
ton amplitude is due entirely to the oscillation which 
took place between the entrance of the Bq region and 
the first bounce. Up to a constant phase factor, it is 
Csi^{asS + apPn). The bounce from the chamber wall 
rotates the amplitude to the new basis (S*, P2) and mul- 
tiplies the individual components by the reflection fac- 
tors As and Ap. Immediately after the bounce, the 



amplitude is -0^"^ = Csi,e^^^{asAsS + apApPn^+i). 
As the particle approaches bounce nL + 1, further os- 
cillation adds CsMflriL-i-i to the amplitude. Immedi- 
ately after bounce nL + 1 the amplitude is ■* = 
Cas{si^e'^^Al + SMAs)S+Cap{si^e'^^Al + SMAp)Pn^+2. 
Repeating this procedure for each of the ng = n^ — nL + 1 
bounces inside the Bq region, we find 



Cas SLe'«M^- + J] U 



n=l 

n_B -1 



+Cap SLG^^Mp" + -SM P^^+^- (50) 



Finally, we add the contribution Csne^^^cin^+i due to 
oscillation between bounce and the exit of the Bq 
region, leading to 



Cas ( SLe'«M"« 



As - 



l-As 
VCap fsLe'«MJ« +SM^ 



Ap - 



V'. 



\B exit) 

7 



A, 



As - A^- 
1- A'; 



(51) 

SRe'«-^ S 
+ SRe*«'^^ P„„+i 
(52) 



srg 



where (5 — > P) denotes the preceding terms with S re- 
placed by P. 

After exiting the magnetic field region, the particle 
bounces from the walls N — nR more times. This mul- 



tiplies the S and P components of ip'] 



Ap~"^, respectively. No further oscillation occurs. Once 
the particle reaches the end of the chamber, the exit 
window performs a quantum-mechanical measurement of 
particle content. Since photons pass through the window 
while chameleons reflect, the chameleon-photon super- 
position is collapsed into one of those two states. The 
probability that this measurement results in a photon is 
"0* • Tpj evaluated at the exit window. Assuming that 
a photon produced by a chameleon with initial condi- 
tions {9, ip) will reach the detector outside the cham- 
ber as an afterglow signal, the expected number of after- 
glow photons generated by this particle is '0* • ipj in the 
small-mixing-angle limit. The time taken by the particle 
to reach the exit window from the entrance is £tot sec 9, 
where £tot = £1 + L + £2 is the total chamber length. The 
contribution of this particle to the detectable afterglow 
rate Faft is the expected photon number over the time, 

I (0f • ^)A^-"-|2cos(0)/4ot + {S^P). 

Only a small fraction of the photons produced emerge 
from the chamber and reach the detector. Most are either 
absorbed in the walls or exit the chamber but miss the 



by 



and 
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detector. In order to find the total decay rate Fdec of 
the population of chameleon particles in the chamber, we 
must account for all of these photons. Since oscillation 
stops after the Bq region exit, and since we already have 
I'ip-yl'^ at that point, we need only compute the probability 
of photon absorption inside the Bq region. At bounce n, 
the absorption probability is given by the difference in 

photon probabilities before and after the bounce, vj;^l = 

l^pj"" ■'p — ''P- Computing the photon amplitudes 
as above, we find 



section 1 



lens 





B field region 
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vacuum chamber 
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^^win,1 










Ifwin 








1 ^^in,2 


B field region 






lens 



I 
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shutter 



V. 



in) 

abs 



-\As\') 
{S ^ P). 



She 



1-A, 



(53) 



The total absorption probability is the sum over all n 
from riL to riR: 



FIG. 9: (Top) GammeV apparatus and (Bottom) CHASE ap- 
paratus (side views, not to scale). Note that CHASE has no 
section 2 {£3 = 0) and has glass windows dividing the Bq 
region into three partitions. Values for these parameters are 
given in Table [H 



(Sexit) 
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C'al{l~\Asn 
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+ SLSM 
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l~\As\' 
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(54) 



The contribution to Fdoc of this particle with initial {6, 
ip) is found by dividing the total conversion probability 

by the time, (IV^^"'^"^^ + T'^f/'^"^) cos(0)/^tot. 

Finally, we average over angles in order to obtain the 
decay rate Fdoc and the afterglow rate Faft: 



abs 



^;(Bcxit) 



COS! 



(55) 



1 



Taft = — / d'n 

Att In 

+{S^P). 



itot 

(56) 



where il/2 is the half-sphere with < 6 < tt/2 and 
< (fi < 2tt. Here, Pdct{d), the probability that the par- 
ticle will reach the detector, is a geometric factor which 
will be computed in Sec. IV Bl The extra factor of 2 in 
the decay rate is due to the fact that particles travelling 
in the negative z direction, away from the exit window, 
contribute to the decay rate but not to the afterglow rate. 
Examining the integrands, we see that the only quantities 
which depend on (p are Og and a|,; each term in the inte- 
grand contains one of these multiplying a (y9-independent 
quantity. Thus the integrals over ip may easily be per- 



formed using 



1 f^'' O 1 p 

— / Updip ~ - cos 0. 



(57) 
(58) 



B. Geometric factor Pdct(6') 

References [H, [13] describe the geometry of the 
CHASE experiment. Diagrams of both experiments are 
shown in Fig. [S] with numerical values given in Table H] 
Note that, in CHASE, £3 = 0; there is no second cham- 
ber section. Henceforth we set £3 = and Ri = R2 = R. 
Outside of the oscillation chamber, a distance ^4 beyond 
the exit window, is a lens with radius ricns and focal 
length ^5. At the focal point of the lens is positioned 
a photomultiplier tube (PMT) whose accepting area has 
radius rpMT- 

We wiU see in Sec. EIIl that the CHASE Bq region 
is divided by glass windows into three partitions. This 
smoothes out sharp features in the oscillation rate during 
the production phase of the experiment. However, we 
will postpone a discussion of these partitions until that 
section, neglecting them for now. 

The previous discussion assumed that the initial 
chameleon particle started out at the center of the en- 
trance window. In the limit that the oscillation length 
AiTio/rn^^ is small compared with the cylinder radius, or 
™cff ^ 0.01 eV in GammeV, the dominant errors in- 
troduced by this approximation are purely geometrical. 
Thus we only need to compute them once, independent 
of the chameleon model. 

The simplest geometric effect arises from the fact that 
the lens and detector are centered on the z axis, which 
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quantity 


GammeV 


CHASE 




2.36 m 


1.61 m 


L 


6.0 m 


6.0 m 




1.16 m 


1.74 m 


h 


2.51 m 


N/A 




2.03 m 


0.64 m 




0.10 m 


0.10 m 


•^win, 1 


N/A 


2.61 


■^win,2 


N/A 


2.91 


Rx 


2.38 cm 


3.175 cm 


R2 


1.75 cm 


N/A 


^win 


N/A 


1.27 cm 




2.54 cm 


2.54 cm 


rpMT 


2.5 mm 


2.5 mm 


to 


1006 sec 


120 sec 


At 


3616 sec 


600 - 3600 sec 


fvci 


0.53 


0.53 


k ^ (jj 


2.33 eV 


2.33 eV 


B 


5.0 Tesla 


5.0 Tesla 


Vpump 


0.026 m^ 


0.014 m^ 



TABLE I: Properties of the GammeV and CHASE experi- 
ments, including the dimensions shown in Fig. [S] 



is where particles start out in our approximation. This 
artificially enhances the probability that a regenerated 
photon will reach the detector. 

A more subtle geometric effect gives particles in our 
approximation a different radial probability distribution 
from the average particle. The volume of cylindrical 
shell at a distance r — {x^ + y^Y^"^ from the z axis, 
AV = 27r£tot''Ar, increases with r. Thus a homoge- 
neous, isotropic gas of chameleon particles will have a 
radial probability distribution P(r) cx r. However, con- 
sider a particle beginning at the center of the entrance 
window; without loss of generality, assume > and 
</j = 0. At the first bounce, and all subsequent odd- 
numbered bounces, x = R. At the second bounce, and 
all subsequent even-numbered bounces, x = —R. The 
path taken by this particle, when projected onto the xy 
plane, will bounce back and forth between the same two 
points. Following [U we call these "2-point paths." All 
r values are equally likely along such a path, implying 
that the radial probability distribution P(r) = constant. 
Therefore, compared with the average particle, the par- 
ticles which we consider spend more time at low r and 
arc more likely to reach the detector. Incidentally, in 
the case of GammeV, particles beginning at r = are 
R1/R2 times more likely to enter the second section of 
the chamber because of this effect. 

We can account for both of these geometric effects 
at once by properly normalizing our afterglow rate in- 
tegrand. The fraction of 2-point paths reaching the de- 



tector is found by integrating over solid angles, 
/2-point = -;— I (f^&dct{0) 

4:71 I SI 

= ^J\medeedetio), (59) 

where Qdct is one if the path reaches the detector and zero 
otherwise. For the CHASE geometry we find /2-point = 
1.57 X lO""*. We do a similar calculation for the average 
particle, whose radial position rex and direction (0ox, V'cx) 
at the exit window are appropriate to a homogeneous, 
isotropic distribution: 

1 1 /•« /-f 

favg = —f^-r- / 27rrcxcfrcx / sm^cx^^'ox 
ttR^ Att /o /o 

X / rf'PoxQdct (^oxj ^oxj V'cx). (60) 

Jo 

For the CHASE geometry /avg = 9.84 x 10"^ Thus, 
when we restrict our calculation to chameleon particles 
beginning at the center of the exit window, we increase 
our afterglow rate by a factor of /2-point//avg = 1.59. 
We correct for this by normalizing P^etiO) in ([55)) : 

Pdctio) = j^^edct((?). (61) 

J2 — point 

In the case of GammeV, we include an extra factor of 
i?2 / ^1 to account for the increased probability of reach- 
ing the second chamber section. This simple, model- 
independent normalization turns out to be reasonably 
accurate. The properly normalized afterglow rate cal- 
culated for 2-point paths differs by 6% from that calcu- 
lated for the "3-point paths" of (SJ] and by 8% from a 
Monte Carlo calculation with arbitrary initial conditions. 
This applies at low masses TOoff ^ (47ra;/L)^/^, or about 
0.001 eV in GammeV. 

Unfortunately there is no equivalent procedure for cor- 
recting the decay rate. We will see that the 2-point de- 
cay rate disagrees with the Monte Carlo calculation of 
Sec. IVIIll by about 40%. However, the experimental up- 
per bounds on the photon coupling P-y are extremely in- 
sensitive to the decay rate; at low /3j, the decay time 
rjg^j. is much larger than the duration of the experiment. 
Furthermore, to the extent that Fdec matters, using the 
2-point calculation is a conservative approximation; (j55p 
computed using 2-point paths overestimates the decay 
rate, hence underpredicts the signal. 

Both the decay rate ([55)) and the afterglow rate ([55)) 
depend on /3y and Bq only through the prcfactor C^. 
Thus Fdcc, Taft oc The rates need only be com- 

puted at one value of the photon coupling and magnetic 
field; a simple rescaling can be used to find the rates at 
other values of these two parameters. Figure [TU| plots 
these rates function of chameleon mass in the 2- 

point computation. Both rates are approximately inde- 
pendent of mass when nies <C 10~^ eV and scale as ni~g 
at nicff ^ lO"'^ eV, consistent with the mass scaling of 
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1e-06 - 



1e-08 - 




0.0001 0.001 0.01 0.1 

effective cfiameleon mass m^^ [eV] 

FIG. 10: Two-point afterglow and decay rates for Bo = 
5 Tesla, 13 j = 10^^ and = 0, for the CHASE afterglow 
experiment (unpartitioned). 



VI. THROUGH THE MAGNETIC FIELD 
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012345678 
distance from entrance [m] 



FIG. 11: Photon amplitude vs. z for a magnetic field 
Ba = 5 Tesla with falloff length Azf = 5 cm. A coupling 
13-y — 10^^ is assumed. Arrows at the right show expected os- 
cillation probabilities for a tophat magnetic field. Low-mass 
chameleons effectively see a tophat field while chameleons 
with TTicff > 0.01 eV see a nearly adiabatic evolution of the 
field. 



A. The real magnetic field 

Thus far we have treated the magnetic field as a 
"tophat" function suddenly increasing from zero to its 
peak value and then dropping again. We have assumed 
a magnetic field -Btophat(-z) = BoQ{z - ii)Q{£i + L - z), 
where &{z) is the step function. The real GammcV and 
CHASE magnetic field falls off from its peak value over 
a distance Azf « 5 cm in a way that is well approxi- 
mated by Btanh(z) = BoTAzi{z-h)TAzi{h+L-z) with 
TAzii^) = (l+tanh(z/Azf + l))/2. For models with oscil- 
lation length ^osc = 47rfc/TO^fj S> Azf, or m^s ^ 0.01 eV 
in GammeV and CHASE, this difference is negligible; the 
field is effectively a tophat function. 

Larger-mass chameleons, however, will see a slowly- 
varying magnetic field in these falloff regions. Their 
oscillation probabilities will differ from tophat function 
approximations. In the extreme case that £osc Azf 
the spatial variation of the magnetic field will represent 
an adiabatic perturbation. A particle which begins in a 
chameleon state and experiences an adiabatically evolv- 
ing magnetic field will return to a pure chameleon state 
once the magnetic field evolves back to zero. Thus = 
in the adiabatic limit. 

Since our integrated solution (|13p is not valid for a 
varying field Btanh, the preceding equation of motion 
must be integrated numerically. Figure [TT] shows the 
results for several chameleon masses, assuming a parti- 
cle trajectory along the z axis. For comparison, arrows 
to the right of the plot show the final probabilities ex- 
pected from the tophat approximation Note that 
the tophat approximation is excellent for masses up to 
the dark energy scale of w 2 x 10~^ cV, and is valid up 
to a factor of order unity at twice that scale. However, for 
rrieff ^ 0.01 eV, Bo{z) rises and falls slowly enough to rep- 
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FIG. 12: Oscillation probability vs. chameleon mass for the 
real magnetic field -Btanh(-z) and the tophat approximation. 
P-y = 10^^ is assumed. 



resent an adiabatic evolution of the background, and the 
oscillation probability is less than the tophat prediction 
by some four orders of magnitude. Figure [T^] shows the 
oscillation probability for Bta,nh{z) and i3tophat(^) over 
a range of chameleon masses. As expected, chameleons 
with TOoff ^ 0.01 eV see a nearly adiabatically varying 
magnetic field. 



B. Window in magnetic field region 

Adiabatic suppression of oscillation can be substan- 
tially reduced by causing a sudden change in the back- 
ground chameleon field inside the Bq region. This is 
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FIG. 13: Photon amplitude vs. z for a magnetic field Bo = 
5 Tesla with falloff length Azf = 5 cm, assuming photon cou- 
pling fi^ = 10^^. Vertical dotted lines show the locations of 
the window and the beginning of the Bq, region. Arrows at 
the right show expected oscillation probabilities for a tophat 
magnetic field, assuming that a measurement is made at the 
window. ( Top) Oscillation probability for a particle beginning 
at the entrance window and measured at the interior window. 
{Bottom) Oscillation probability for a second particle begin- 
ning on the interior window is added to the probability due 
to the first particle. 



accomplished by placing a glass window inside that re- 
gion. Close to the surface of this window, the background 
chameleon field and mass (|5T|) change on distance 
scales of order T?T,off(0buik(Pmat))~^- By the containment 
condition this must be less than the wavelength ^ k~^, 
which is in turn much smaller than £osc- 

Since an increase in the chameleon mass sharply 
decreases the chamclcon-photon oscillation probabil- 
ity p^ . the two are effectively decoupled near the 
window. Thus the sudden change in the background 
chameleon field due to the window does not lead to 
any sudden changes in the oscillation amplitude. Fig- 
ure [13] ( Top) shows the numerically computed amplitudes 
for a range of masses up to the location of a window 
placed one meter inside the Bq region. Arrows at the 
right show the corresponding probabilities for a tophat 



magnetic field. The tophat approximation is clearly bet- 
ter here than in Fig. [T^l Of course, the smooth tran- 
sition at the entrance to the Bq region still has some 
effect on transition probabilities. Since it adds a few ex- 
tra centimeters to the length of the Bq region, it shifts 
the locations of the features in the afterglow rate seen in 
Fig.[TUl Also, some of the suppression due to the smooth 
transition remains. This is best seen in the plot for the 
TOoff = 5 X 10""^ eV chameleon in Fig. [13] (Top), which 
has a final photon amplitude about 40% less than the 
tophat prediction. 

Section [V] began with a particle in a pure chameleon 
state immediately after measurement by the entrance 
window of the vacuum chamber. In the case of a win- 
dow in the interior of the Bq region, a particle can also 
begin in a pure chameleon state on the surface of the win- 
dow. By the homogeneity and isotropy of the chameleon 
population we expect that for every particle incident on 
the window from low z with momentum k there will be 
another leaving the window on the other side with the 
same momentum. In order to compute the total oscilla- 
tion probability due to both partitions of the magnetic 
field, we must also consider the oscillation due to this 
second particle. Fig. [T3] {Bottom) shows this combined 
probability, which once again agrees well (with errors at 
the ^ 10% level at high masses) with the tophat magnetic 
field calculation. 



C. Interior windows in CHASE 

CHASE uses two interior windows, the first 1.0 m in- 
side the Bq region and the second after another 30 cm. 
Each window is 1 inch in diameter. Figure [HI f Toy] com- 
pares the resulting oscillation probability to its tophat 
approximation. The adiabatic suppression of oscillation 
by many orders of magnitude seen in Fig. [T^] is not in 
evidence here. However, some suppression remains, as 
shown in Fig. [T3] {Bottom). The points in that figure 
use probabilities averaged in bins of width 10% in rrics 
in order to average out features due to different effec- 
tive Bq region lengths. This averaging shows that the 
oscillation probability for the tanh field is 30% lower 
than the tophat approximation at high masses. Since 

1/2 

cx V^[^^, this implies that constraints on found 
using the tophat approximation will overestimate actual 
constraints by w 15% for rricfr > 0.01 eV. 

A more accurate numerical computation would also 
consider: 

• a range of particle directions fc, since a particle not 
propagating parallel to the z axis would see an ef- 
fective magnetic field falloff length larger than Azf; 



• a range of particle frequencies within the 
lincwidth of the laser; 



1 cm 



• reflections from the chamber walls inside the Bq 
region, leading to partial measurements of particle 
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FIG. 14: (Top) Oscillation probability vs. chameleon mass 
for tophat and tanh magnetic fields, assuming two interior 
windows and /3.y — 10^^. (Bottom) Ratio of tanh and tophat 
oscillation probabilities. 



content. 

Each of these corrections will serve to smooth out fea- 
tures in the oscillation probabilities in Fig. [TH Since nu- 
merically computing the probability for each path used 
in a calculation such as that in Section |V] would be pro- 
hibitive, and since adiabatic suppression is only impor- 
tant in the least interesting region rricff > 0.01 eV of our 
parameter space, we have retained the tophat magnetic 
field approximation in our computations of decay and af- 
terglow rates. Thus we expect that our constraints will 
overestimate actual constraints on /3-y at the ten percent 
level, or approximately the thickness of the line used to 
plot constraints in Figure [25] of Section IIX Dl 



VII. ENHANCEMENTS IN CHASE 

The GammeV chameleon experiment was similar to 
the idealized afterglow experiment discussed in Sec. Ill El 
CHASE made major improvements to the vacuum sys- 
tem and the detector, as well as adding data runs with 



lower magnetic field values and partitioning the magnetic 
field region into sections. In this section we discuss each 
of these improvements. 



A. Vacuum system 

The range of chameleon potentials probed depends on 
the density ratio between the vacuum and the walls of the 
vacuum chamber. A chameleon particle will be trapped 
inside the vacuum chamber if the containment condition 
(PH]) is satisfied for E = 2.33 eV and pmat equal to the 
lowest density in the chamber walls. Meanwhile, an os- 
cillation experiment will only be sensitive to chameleons 
below some mass TOmax at the density pvac of the cham- 
ber vacuum, since massive chameleons mix very weakly 
with photons. 

Here we study power law and chameleon dark energy 
potentials (|8l9p in order to find the range of exponents 
n for which some region of the (/3-y, g) parameter 
space is accessible to a given afterglow experiment. Note 
that our goal in this section is not to compute actual 
constraints, but only to provide a simple estimate of the 
range of potentials which may be probed. It will be useful 
for this purpose to neglect pc,EM, since a greater total 
density inside the chamber degrades constraints. This 
assumption is equivalent to <g; /^m- 

For potentials (|8l9p . assuming negligible pc,EM, 
eq. (|lip for the chameleon mass implies that nics oc 

Pmat^^^'^""^'- Requiring that mcff((?!)buik(pmat)) > E 
(containment) and (^!)buik(Pvac)) < m^a.^ (mixing), 
we find 



2 ^ l0g(£;/TO,nax) 



2n-2 log(pniat//Ovac) 



(62) 



Now we must consider n > 2 and n < chameleons 
separately. If n > 2 and pmat//Ovac < (-B/m-max)^, then 
the containment and mixing conditions cannot both be 
satisfied; the experiment is insensitive to these models. 
If n > 2 and Pmat/Pvac > (£^/"imax)^ then 



n > 2 



l0g(Pmat/Pvac) - log(-E/mniax) 
l0g(/9mat/pvac) " 2 log(iJ/TOinax) 



(63) 



That is, for n satisfying this condition, the experiment 
will be able to probe some region of the (/3-y, /3^, g) pa- 
rameter space. If n < and -E/mmax < Pmat/Pvac < 
(£;/mmax)^ then 



n > 



, l0g(Pmat/Pvac) - log(£^/mmax) 
2 log(i?/TOmax) - log(pmat/Pvac) ' 



(64) 



If Pmat/Pvac > (-B/m-max)^ then all n < can be probed 
by the experiment. 

GammeV used a turbomolecular pump in combination 
with a roughing pump. Since any chameleon particle 
entering the roughing pump would be scooped out of 
the chamber, GammeV was sensitive only to chameleons 
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massive enough to bounce from the P = 1.9 x 10~^ torr 
gas at the intake of that pump. Thus the lowest- 
density "wah" in GammeV had pmat = 4 x 10~^ g/cm^ . 
The vacuum density was pvac = 2 x 10~^^ g/cm^ and 
TOmax ~ 10~^ eV [31[. Thus GammeV was completely 
insensitive to potentials with n > 2 and could probe only 
-0.77 < n < 0. 

CHASE maintained a vacuum of pvac — 3.7 x 
10~i^ g/cm^ using ion pumps and cryo-pumping on the 
4 Kelvin walls of the magnet bore. This vacuum system 
did not allow chameleons to escape from the chamber, 
meaning that the lowest-density regions of the chamber 
wall were the entrance and exit windows, made of BK7 
glass, with density pmat = 2.51 g/cm.^^. Although the ex- 
perimental sensitivity to high-mass chameleons depended 
on the photon coupling, TOmax ^ 0.01 eV is a reasonable 
approximation. With this approximation, CHASE can 
probe 77, < and n > 2.5. However, we note that 
Pmat ^ Pc,EM in CHASE, so the range of potentials 
probed will be smaller when (3-y > /3m. For example, 
when = lO'^Pmat/pc.EM, CHASE is sensitive to 

n < and n > 2.8. 



B. Multiple magnetic fields 

In order to extend constraints on the photon coupling 
to lower values, an afterglow experiment would need to 
use a large magnetic field and collect data over a long 
time period, as was done by GammeV. However, the 
afterglow signal ([T7)) falls off exponentially with a rate 
Fdec oc P^Bq. Thus increasing Bq makes the experiment 
insensitive to chameleon models with large photon cou- 
plings. Since the signal falls off rapidly, collecting data 
for a longer time does not help. 

CHASE used short runs with low magnetic fields to 
supplement the high-_Bo, long-duration runs. For each 
photon polarization, data were collected in two long runs 
of Bo = 5.0 Tesla and filling times of tprod = 5 hours. 
Short data runs of tprod = 10 minutes each were con- 
ducted for magnetic fields of Bq = 2.2 T, 1.0 T, 0.45 T, 
0.20 T, 0.089 T, and 0.050 T. Additionally, seven calibra- 
tion runs were carried out with tprod = 10 minutes and 
i?o = to allow for further study of backgrounds. 

These low-i?o runs extended constraints to models with 
larger . Figure [15] shows the afterglow signal Faft (t) 
from (|17[) for various magnetic fields and photon cou- 
plings. The narrower, light blue shaded region represents 
the time window used for the short runs. The larger, yel- 
low shaded region is the extra observation time in the 
long runs. A detailed analysis of the data, properly ac- 
counting for statistical and systematic uncertainties, will 
be provided in Sec. IIXI Here, Fig. [T5] roughly approxi- 
mates CHASE constraints by assuming that an afterglow 
signal of at least 1 Hz two minutes after the laser is turned 
off will be detected by the experiment. Thus a curve in 
Fig. [15] entering the shaded regions will be detected. 

At the lowest fi^ in Fig. [THl the chameleon is detectable 



only in the 5 Tesla run. By (3^ = 10"'^^ several other runs 
can also probe the chameleon. Meanwhile, the decay time 
T'l^ in the 5 Tesla run is less than the observation time, 
so that a decline in i^aft is apparent. At = 10^^ the 
signal in the 5 Tesla run falls off so rapidly that it is no 
longer observable, and CHASE constraints are due en- 
tirely to the short runs. By (3^ = 10^^, Fdec is so large, 
even in the low-_Bo runs, that none of the runs can de- 
tect the chameleon signal. Since collider constraints [i^ 
already exclude > Mpi/(lTeV) = 2 x 10^^, further 
CHASE runs at still lower magnetic fields are unneces- 
sary. 

C. Interior windows 

As described in Sec. IVI C[ the magnetic field region 
in CHASE is divided into three partitions by interior 
windows. In addition to reducing the adiabatic suppres- 
sion of oscillation rates, these windows fill in some of the 
gaps in the constraints due to destructive interference 
during the production phase. Without the partitions, 
the afterglow signal (fTT)) contains a factor of V-^^ip due 
to chameleon production as photons stream through the 
chamber. V^^^cf, vanishes when m^^ = Anirk/L for any 
integer n. This corresponds to maximal destructive inter- 
ference in oscillation; the chameleon amplitude produced 
in one part of the chamber exactly cancels that from an- 
other. 

Figure [16] shows the chameleon production probabil- 
ity per photon, 7'-y<^0, as a function of effective mass, 
with and without partitions in the magnetic field re- 
gion. Without partitions, the probability has several ze- 
ros near the dark energy mass scale of 2.4 x 10^"^ eV. 
The partitions push the first zero of 7^-y<-i.0 from rrtcs ~ 
9.8 X 10""* cV to 7.6 X 10~'^ eV, eliminating the sharp 
features in the predicted signal in the region of the dark 
energy scale. In fact, suppressing destructive interference 
was the original reason for inclusion of interior windows 
in the chamber. The benefits of these windows for mit- 
igating adiabatic suppression of chameleon production 
were discovered only later. 



D. Detector 

The dominant source of systematic uncertainty in 
GammeV was instability in the photomultiplier tube used 
to detect afterglow photons [3l|. The dark rate seen 
in the PMT had large fluctuations on time scales of 
1 minute. As a result, constraints on chameleons had 
to be obtained by averaging the signal reported by the 
PMT over long times. Even the time- averaged rate had 
a substantial uncertainty; averaging over a large num- 
ber of hour-long intervals resulted in a mean dark rate of 
115 Hz with a standard deviation of 12 Hz. 

CHASE used a PMT with a lower mean dark rate, 
30 Hz, which was much more stable with time. More- 
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time [sec] time [sec] time [sec] time [sec] 

FIG. 15: Afterglow signal vs. time for various magnetic fields (labelled in plots) and photon couplings (upper right corners of 
plots). Approximate observation windows for the short data runs (thin, light blue shaded regions) and long data runs (both 
shaded regions) are shown in each plot; the long runs are used only for Bq — 5 Tesla. In all cases, rries ~ 10~* eV is assumed. 
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FIG. 16: Production probability for a chameleon with I3-, = 
10^^ in a magnetic field Bo = 5 Tesla, with and without the 
CHASE partitions. 



over, the PMT was modulated by a shutter with a period 
of 30 sec and a duty cycle of 0.5, allowing minute-scale 
variations in the dark rate to be monitored in real time. 
Fortunately the dark rate of the CHASE PMT remained 
stable. Nevertheless, the experiment was designed to de- 
tect any such instability and to subtract its effects from 
the signal. 



VIII. CHAMELEON-PHOTON OSCILLATION: 
MONTE CARLO SIMULATION 

Section |V] computed the afterglow rate Faft and de- 
cay rate Fdec per chameleon particle. That calculation, 
which assumed a simplified set of initial conditions, was 
fast and accurate for a GammeV-like experiment. Its 
extension to the CHASE geometry was found to be ac- 
curate to a factor of ~ 2 over most of the mass range, 
makeing it useful when designing broad features of the 
CHASE geometry. However, it was inadequate for data 
analysis for several reasons. 

1. The assumption that the chameleon particle's ini- 
tial position on the entrance window is irrelevant 



breaks down when the oscillation length 47rfc/m^fj 
of the theory becomes smaller than the chamber 
radius, corresponding to a mass of rrics ~ 0.01 eV. 

2. Even at low chameleon masses, the calculation of 
Sec.|V]is inaccurate for nontrivial chamber geome- 
tries such as the interior windows in CHASE. 

3. That calculation cannot include random systematic 
effects such as diffuse reflection from the chamber 
walls. 

Thus a Monte Carlo simulation of CHASE was used to 
compute the afterglow signal and to estimate the effects 
of random systematics in the CHASE analysis [s^] ■ 



A. Particle paths 

Relaxing the restrictions on the initial conditions used 
in Sec. |Vl the decay and afterglow rates are given by 
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where the (Px integrals are carried out over the surfaces 
of the windows and the (PQ integrals over the unit sphere 
of k values. Here Awin is the total surface area of the win- 
dow surfaces inside the chamber. The probabilities Paft 
and Pdcc are oscillation probabilities along particle paths 
which are completely determined by their initial positions 
and directions. A Monte Carlo calculation of these rates 
replaces the integrals with sums over particle paths with 
randomly chosen initial conditions. In the limit that the 
number of particles simulated becomes large, the Monte 
Carlo rates will approach the actual rates. 

A particle in a pure chameleon state can result from 
any quantum-mechanical measurement of particle con- 
tent at a glass window, including the entrance and exit 
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windows as well as the interior windows of CHASE. Thus 
we initialize each particle by randomly choosing a lo- 
cation on one of the window surfaces, with probability 
proportional to surface area. When computing the de- 
cay rate, the initial particle direction is chosen randomly 
from a uniform distribution as appropriate to an isotropic 
chameleon population. Although this would also be cor- 
rect for the afterglow calculation, photons emerging from 
the chamber at most of these directions would not reach 
the detector. This would cause the Monte Carlo calcu- 
lation of Faft to converge very slowly. Instead we choose 
particle directions from a distribution that is uniform 
over the portion of the unit sphere with 6 < 9o = 0.1 
and zero elsewhere. Since this probability distribution 
covers a fraction = (1 — cos0o)/2 of the unit sphere, 
we must multiply our result by fe to obtain Faft- 

Once an initial position and direction have been cho- 
sen, the particle is propagated in that direction until it 
encounters a wall or window. At a chamber wall, the 
particle direction is changed and its photon momentum 
corrected to account for absorption and phase shifting, 
as described in Sec. IIVI The particle is then propagated 
forward once again until a window or another wall is 
encountered. At a window, the particle content is mea- 
sured. Continuing to follow the chameleon particle af- 
ter this point would constitute a double-counting, since 
the chameleon position and direction immediately after 
this bounce from the window is within the allowed set 
of chameleon initial conditions. For this chameleon path, 
Pdoc is the probability that the particle is a photon at the 
window plus the total probability of absorption during 
collisions with chamber walls. In order to compute Paft, 
we must compute not just the photon amplitude at this 
window but the probability that the photon will escape 
from the chamber and reach the detector. The factor 
|^^|2(Ar-nR) ^ probability of escape without 

absorption in the walls, and the factor Pdet is the proba- 
bility of detection. Thus we must also keep track of the 
photon which may result from quantum measurement at 
a window. Figure [T7] shows a few sample paths used by 
the Monte Carlo calculation to compute the afterglow 
rate (Left) and the decay rate (Right). 

From Fig. [17] it is evident that simulated paths used 
in the decay rate computation bounce many more times 
from the chamber walls than paths used to compute the 
afterglow rate. This is because paths which bounce many 
times have large angles 9 with respect to the cylinder axis; 
thus photons resulting from such paths are less likely to 
emerge from the CHASE chamber and to reach the de- 
tector. Figure [TS] shows the fractional contributions to 
the total afterglow and decay rates from paths with dif- 
ferent numbers of wall bounces. As expected, afterglow is 
dominated by paths with very few bounces; a third of the 
paths bounce exactly twice, while over 90% bounce be- 
tween 1 and 5 times. A path with only two bounces must 
have an ang le less than 2i?/(i4ot) = 0.020, or about 1°. 
Thus the afterglow rate is dominated by paths reflecting 
from the chamber walls at grazing incidence. In con- 
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FIG. 17: Typical particle paths used in the Monte Carlo sim- 
ulation of CHASE. (Top) Paths used to compute Faft. The 
red path (with "-f " symbols showing bounces from chamber 
walls) begins on an interior window and ends on the exit win- 
dow. The green path (with "x" symbols denoting bounces) 
begins on the entrance window and ends on the exit window. 
The blue path (with "*" symbols denoting bounces) begins 
on the entrance window and ends on an interior window; the 
subsequent dotted line denotes the photon which may result 
from the quantum measurement at the interior window. {Bot- 
tom) Paths used to compute Fdec- The red path ("-I-") begins 
on an internal window and ends on the entrance window, while 
the green path ("x") begins on the exit window and ends on 
an interior window. Since a larger range of directions is al- 
lowd in the decay rate computation, the typical path bounces 
from chamber walls much more frequently than those used in 
the afterglow computation. 



trast, paths contributing to the decay rate reflect at a 
large range of angles. Their contributions to the total 
rate do not fall off until the number of bounces is around 
£tot/i? « 300. 



B. Tests of the Monte Carlo simulation 

In the limit that the number iVparticies of particles sim- 
ulated becomes large, the Monte Carlo computations of 
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FIG. 18: Contributions to (Top) afterglow and to (Bot- 
tom) decay rates from paths vs. number of wall bounces per 
path. The afterglow rate receives no contributions from paths 
bouncing more than 9 times, since these emerge from the 
chamber at too great an angle to reach the CHASE detector. 
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FIG. 19: Convergence of the Monte Carlo simulated rates 
as the number of particles simulated is increased. For each 
particle number A''particies , each point represents a different 
random number seed used in the Monte Carlo simulation. 
= 10^^, mcff — 10^* eV, and Bo — 5 Tesla are assumed. 



tides begin at the center of the entrance window. The 
same conditions can be imposed on the Monte Carlo cal- 
culation for the sake of comparison. Fig. [20] shows Taft 
and Tdoc computed using the two methods. Evidently 
both rates are accurate at the 1% level across the full 
range of chameleon masses. 

Figure [H] shows the afterglow and decay rates per 
chameleon particle, computed with no restrictions on the 
initial particle position, and including the interior win- 
dows. In Section llXl wc will use such Monte Carlo calcula- 
tions with appropriate values of the potential-dependent 
phase £,v in order to determine CHASE constraints. 



Taft and Fdec should converge to constant values. This 
convergence can be studied by changing the random num- 
ber seed used by the simulation at fixed -/Vparticics- The 
scatter in the resulting rates should be smaller at larger 
-^particles- Figure fT9l shows that the rates do indeed con- 
verge, in the sense that the standard deviation of several 
rate calculations with different random number seeds be- 
comes small. The decay rate Fdcc has converged to 0.5% 
by iVparticies = 10^ and the afterglow rate has converged 
to 1% by A^particies = 10^. Sincc most of the simulated 
afterglow paths do not reach the detector, even with the 
restricted set of initial conditions used, -/Vparticics must be 
an order of magnitude larger in the afterglow calculation 
in order for the rates to converge to the percent level. 
Henceforth we use iVparticios — 10^ for the Fdcc computa- 
tion and -/Vparticics = 10^ for the Faft computation. 

We have shown that the Monte Carlo computation of 
Faft and Fdoc is precise. Next we show that it is accu- 
rate. The 2-point calculation of Sec. |V] found these rates 
exactly for a chamber with no interior windows and a 
restricted set of initial conditions, namely, that the par- 



C. Diffuse reflection 

Thus far we have assumed perfectly specular reflection 
of particles from the CHASE chamber walls; that is, the 
angle of incidence equals the angle of reflection. Sup- 
pose instead that a fraction f^is of bounces resulted in 
perfectly diffuse reflection, in which all directions were 
equally likely for the reflected particle. This could result 
from a rough surface whose local normal vector could dif- 
fer substantially from that of a perfect cylinder. Given 
an incident direction kj and a randomly chosen reflected 
direction kn, the effective local normal vector is propor- 
tional to kfj^ kj. 

Figure [22] shows the change in afterglow and decay 
rates when the fraction of diffuse reflections is increased. 
For this computation we allow all initial directions for 
afterglow paths, rather than requiring that 6* < 6'o = 0.1 
as before. The qualitative effect at fdis ^ 1 is that dif- 
fuse reflection increases the afterglow rate and slightly 
decreases the decay rate. This is because such diffuse 
reflection is most important for particles with large 9, 
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FIG. 20: (Top) Taft and Tdcc for the 2-pomt (2pt) calcula- 
tion of Sec. |V] as well as the Monte Carlo (MC) calculation 
with all particles required to start at the center of the en- 
trance window. In both calculations, the chamber has no 
interior windows. /?-y — 10^^ and Bo = 5 Tesla are assumed. 
{Bottom) Fractional difference between the 2-point and Monte 
Carlo calculations. 



which bounce many times but do not reach the detec- 
tor. Diffuse reflection gives such a particle a nonzero 
chance of reaching the detector, thereby increasing the 
afterglow rate. Since increasing Faft and decreasing Fdec 
both improve CHASE constraints, henceforth we make 
the conservative assumption that fdm = 0. 



D. Sensitivity to chamber properties 

The geometry of the CHASE experiment is illustrated 
in Figure [9l with numerical quantities listed in Table ID 
The length L of the magnetic field region, as well as its 
offsets £i and £2 from the entrance and exit windows, 
respectively, affect the afterglow and decay rates. The 
same is true of the reflectivity of the chamber walls. Such 
properties may be important to the design of future af- 
terglow experiments, so it is instructive to consider their 
effects on the CHASE rates. 
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FIG. 21: Afterglow and decay rates computed using the 
Monte Carlo simulation. Particles are allowed to begin at 
any position on a window surface, = 0, = 10^^, and 
Bo = 5 Tesla are assumed. 
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FIG. 22: Fractional change in afterglow and decay rates 
vs. fraction fdis of diffuse reflection. Rates for six different 
random number seeds are shown. 



Figure ( Toy) shows how Faft and Fdcc change when 
the lengths £1, L, and £2 are varied relative to their 
CHASE values. As expected, increasing L from low val- 
ues causes Faft and Fdec to grow as L^. Particularly for 
the afterglow rate, this growth levels off as destructive 
interference during oscillation becomes important and as 
each particle must bounce from the walls a greater num- 
ber of times before reaching the exit window. Mean- 
while, increasing £1 and £2 cause the chamber volume to 
increase without adding volume to the magnetic field re- 
gion. This means that each particle spends less time in 
the Bq region, hence the oscillation rate is smaller. In- 
creasing £2 also has the effect of making particles bounce 
more between the Bq region and the exit window, further 
suppressing Faft- However, varying £1 and £2 by factors 
of two in either direction only changes Fdcc by about 10% 
and Faft by at most 25%. Thus CHASE constraints will 
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FIG. 23: (Top) Afterglow rates (thin lines) and decay rates 
(thick lines) as the lengths £i, L, and £2 are varied relative 
to their CHASE values in Table |I] {Bottom) Afterglow and 
decay rates vs. the mean wall reflectivity frat in the thin-skin 
limit hi —>■ 00. In both plots, meff — 10~* eV and ~ 
have been assumed. 



be relatively insensitive to the chamber geometry. 

An additional consideration is the volume occupied 
by the vacuum system. Suppose that the cylindri- 
cal chamber considered above is only a fraction /voi of 
the total volume available to chameleon particles. For 
CHASE /vol = 0.68. Since the distribution of chameleon 
particles is homogeneous, only a fraction /vol of the par- 
ticles will be in the chamber at any given time, with the 
remainder inside the vacuum system. Thus the rates Tdoc 
and Faft computed earlier in this section must be mul- 
tiplied by /vol. Constraints presented in Sec. IIXI include 
this factor. 

The mean reflectivity /rof is varied in Figure [23l (Bot- 
tom). Since the photon skin depth and the visibility fac- 
tor /vis discussed in Sec. IIV Cl are specific to the material 
of the chamber wall. Fig. [23] {Bottom) makes the simpli- 
fying assumption of zero skin depth. This is equivalent 
to ni — )• oo and /vis = fid in the notation of Sec. IIV CI 
Around the CHASE value of /rcf = 0.53, the afterglow 
and decay rates scale approximately linearly with /rcf. 



We note that the 7.5% measurement uncertainty in /j-of 
will lead to « 10% uncertainties in Faft and Fdoc- Since 
(X i^aft^ ^ uncertainty in Faft implies 

a 2.5% uncertainty in the CHASE upper bound on /3^, 
which is a nontrivial contribution to the total uncertainty. 

In future experiments, highly polished chamber walls 
would not strengthen constraints by much. Polished 
metal, with /rcf ~ 0-85 — 0.9, would only double the af- 
terglow signal, improving constraints on fi^ by 2^/^ — 1 ss 
20%. What is important is measuring /rcf at the 10% 
level so that a precise bound can be placed on 



IX. ANALYSIS AND CONSTRAINTS 
A. Profile likelihood analysis 

Finally, we apply the results of the preceding sec- 
tions to CHASE data using the profile likelihood method 
of [3^ . The afterglow photon rate i^aft from ([TT]) depends 
on the afterglow and decay rates per chameleon, Faft and 
Fdcc, as well as the production time tprod- The rates Faft 
and Fdcc, as functions of the chameleon parameters rrics 
^v, and Pj, for CHASE as described in Sec. IVIII are 
computed in Sec. IVIIIl Fig. [5l] shows these rates for a 
particular choice of and /3^. Thus for CHASE we 
know Faft {rUeS ,£,v,P-i,t). 

A diagram of CHASE is shown in Fig. |9] {Bottom). In 
each run, during the afterglow phase of the experiment, 
CHASE counts photons in 15 second bins. A shutter 
covers the detector, a PMT, in every other bin, allowing 
the background photon "dark rate" to be monitored in 
real time. The excess photon rate observed by the PMT 
must come from the vacuum chamber. It is some com- 
bination of the afterglow rate Fatt(?Tioff , ^y, /3-y, t) and a 
background systematic rate Fgystlipi}, which depends 
on "nuisance parameters" {pi}. 

Let the eight data runs and seven calibration runs be 
labeled by r, and the time bins in run r by 6^; the mean 
photon rate observed in the fo^th bin of run r is F^f^^'^ . 

The predicted photon rate f'J^"'^'(TOcff , , f^-y 1 {pi}) , for 
a given choice of the chameleon and nuisance parame- 
ters, is found by averaging Faft + ^syst over time t in 
the appropriate bin. The uncertainty cr^.h^ in that bin is 
d ominated by the Poisson noise in the 28 Hz dark rate, 
■\/28Hz/15sec = 1.4 Hz, and contains additional noise 
from background photons (32l |. Summing over all runs 
and all bins in each run, we define 

/ ^(obs) _ ^(prcd)\^ 

x\m,s,^v,[3j,{p^}) = J2 '^'''^ 2 '^''^ ■ (69) 

At each point in the chameleon parameter space, the 
profile likelihood method defines x^{'mcji,£,v : P-i) to be 
X^("^ofT, 'Cv J /37J {pi}) minimized over the nuisance pa- 
rameters. This is compared with the value Xnuii fo'" 
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null model = -Faft = 0, which has no photon-coupled 
chameleon field. 



B. Systematic rate Fsyst 

The total background rate Fgyst is found to have three 
important components [32|: 



1. a dark rate Fdark = 28 Hz in the PMT, which is 
the same for all data runs; 

2. a "glow" Fpump emitted by the ion pump, which 



varies from run to run with a mean of Fi 
1.2 Hz and a standard deviation of a^, 



pump, mean 



0.4 Hz; 



3. a transient rate i^tran(i) which is the same function 
of time for all data runs. 

-Fdark IS degenerate with individual pump glow rates, so 
we fix it in the analysis. 

At late times, the background is dominated by the 



first two of these, F=vst —J- Fd 



"syst ^ J^dark ~t~ -^pump ^- kMUCC dark 

are time-independent, neither one depends 
upon the bin number br- Poisson variations in i^dark of 
\/ Fda.i±/ Atr^br the dominant component of the uncer- 
tainty in a bin of width At^^b^. The set {pi} of nuisance 
parameters includes the pump glows Fpu^p with the 
mean and standard deviation above, but not the known 
dark rate i^dark- We modify (|69)) above to include the 
term J^A^p 

uinp ^ -^puinp,incan 

)^/'^pump in order to ac- 
count for the uncertainty in Fpump ^. 

The transient rate Ftran (t) is studied in [111 . It cannot 
be a chameleon afterglow because it is independent of 
the magnetic field, and because its amplitude peaks in 
the orange region of the electromagnetic spectrum rather 
than the green 2n/uj = 532 nm of the input photons. 
Because of this spectrum, references [13, l35| refer to this 
transient component as the "orange glow" . It is modeled 
as a run-independent exponentially decaying background 
photon rate Flran(i) = -For.o exp(-ror,oi)- 

Although i^tran canuot be a chameleon afterglow sig- 
nal, it can mimic such a signal in one particular run. 
Furthermore, -Ftran will introduce correlations among the 
errors in the first several data bins. Therefore we cannot 
simply measure i^tran and subtract it from the observed 
signal. We must treat it as a systematic to be fit; that 
is, we must include i^or,o and For.o in the set {pi} of 
systematics parameters. Analysis of CHASE calibration 
data shows i^or,o ~ 7 Hz and l/ror,o ~ 2 minutes. In or- 
der to avoid the portion of the data most contaminated 
by this systematic, henceforth we discard the first two 
minutes of afterglow data in each run. 

The "profile x^" minimizes over the nuisance param- 
eters {pi}- Since the total systematic rate Fsyst{t) de- 
pends linearly on all of the nuisance parameters except 
for For.Oj given a value for Fqi-.o the x^-minimizing values 
for the others can be found by solving a linear system. 



Since F^i 



and Fr 



pump - 



Letting t^.b^ and Atr.b^, respectively, be the central time 
and the width of bin br in run r, we define 



Jtr,b,.-At^,i,^/2 



(70) 



Sr,br = ^Sinc f iFor,oAir,6, ) CXp(-ror,oir,bJ (71) 
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/p(obs) _ 77'(aft)\<^TO 

n _ \ " ^ r,br r,br- )^r,l 



(72) 
(73) 

(74) 



br 



i-,b,. 



for m = 0, 1, 2. Then for.o and the i^pump ^ are given by 



2 ~ A ' 



(75) 



F 



pump ^ 



2F^or,o-A^. (76) 



•^r,pump 



•^r,pump 



must be minimized numerically with respect to the 
remaining nuisance parameter, For,Oj in order to deter- 
mine the profile x^, which depends only upon moff, ^y, 
and /3-y. Henceforth we will use to refer to the profile 

defined above. 



C. What would a chameleon signal look like? 

Simulated CHASE "data" show what a chameleon 
afterglow signal would look like in the experiment. 
For a given chameleon model, the predicted rate 

Fl^,^'^^\mce,£,v, [3j,{pi}) was simulated by setting 
-For.o = 7 Hz and For,o = 1/120 sec, and by choosing 
the pump glows for each run from a Gaussian distribu- 
tion of mean 1.2 Hz and standard deviation 0.4 Hz. In 
bin br of run r this implies F^'^^'^'^'' Atr^br photons. The 
simulated "data" for this bin are generated by randomly 
choosing a number from a Poisson distribution of mean 

"Data" were simulated for three different scenarios: 

1. a chameleon with trics = 10^* eV, £,v = 0, and 

~ 10^'*, near the center of the CHASE sensitiv- 
ity region; 

2. a chameleon with TOoft = 2.4 x 10"'^ eV, = 0, and 
I3.y = 10^^, near the edge of the sensitivity region; 

3. a null model, with no chameleon afterglow. 

Figure [51] shows the constraints resulting from analysis 
of these three simulations. For the first scenario, ana- 
lyzed in F\g.^^{Left), CHASE data would be extremely 
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FIG. 24: Analysis of simulated CHASE data. Yellow regions show previous constraints; blue regions show models excluded 
relative to the null (no-chameleon) model by Ax^ ~ 6.0, corresponding to a 95% confidence level for a Gaussian probability 
distribution; green regions show chameleon models preferred relative to the null model by Ax^ ~ 9.2, corresponding to a 99% 
confidence level; pink regions within the green regions show chameleons within Ax^ = 6.0 of the best-fit chameleon model. 
Black signs show the fiducial chameleon model. (Left) Simulated chameleon near the center of the CHASE sensitivity 
region, tticB — 10^'' eV and Pj = lO^**. The inset shows 0.99 x 10~^ eV < tticB < 1.02 x lO^** eV on the horizontal axis and 
0.9999 X 10^"* < < 1.0001 x lO" on the vertical axis; the black ellipse shows the Ax^ = 2.3 contour, corresponding to the 
68% confidence level. (Middle) Simulated chameleon near the edge of the CHASE sensitivity region, nies = 2.4 x 10~^ eV and 
— 10^^. (Right) Simulated data with no chameleon. In all of the analyses, = has been assumed. 



powerful. The experiment would be able to constrain the 
chameleon mass to '--^ 1% and the photon coupling to bet- 
ter than 0.01%. In the second scenario, CHASE would 
detect the presence of a chameleon to high significance 
(Ax^ > 50), but severe parameter degeneracies would 
make a determination of nicff and (5^ difficult. Fur- 
ther study of this chameleon would require a new or 
redesigned experiment. Finally, in the third scenario, 
CHASE would exclude chameleon models over a large 
range of parameters. The profile likelihood analysis pre- 
vents a spurious identification of the orange glow system- 
atic as a chameleon afterglow. 



D. CHASE model-independent constraints 

Now that the CHASE analysis has been tested on simu- 
lations, we proceed to the actual data. Here we calculate 
constraints which are model-independent, in the sense 
that the chameleon parameters are the effective mass nias 
inside the oscillation chamber, the phase ^v, and the pho- 
ton coupling /3y , rather than the potential and the matter 
coupling. In order for these constraints to be applicable, 
the chameleon mass in the chamber walls must be large 
enough that chameleon particles are contained inside the 
chamber. 

Figure shows CHASE constraints [s^ on scalar and 
pseudoscalar chameleons with = 0. Yellow regions 
are previous constraints in this parameter space; Gam- 
meV constraints are taken fro m I Sll . and collider con- 
straints were found by [1^ H^, l4l|. In the blue region, 
the null model is preferred relative to the chameleon 
by Ax^ = 6.0, which corresponds to exclusion at the 
95% confidence level for a Gaussian probability density 
function; henceforth we use "95% GL" to refer to this 
Ax^ = 6.0 contour. CHASE constraints improve upon 
those of GammeV by: extending to higher mcff, well be- 
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FIG. 25: Model-independent chameleon constraints for = 
0, from [s^]. The solid blue region is exluded at the Ax^ ~ 6.0 
level (95% CL for Gaussian probability) for scalar chameleons, 
and the interior of the green curve is excluded at that level 
for pseudoscalar chameleons. 



yond the dark energy scale of 2.4 x 10^"^ eV; bridging the 
gap between GammeV and collider constraints; improv- 
ing the low-moff upper bound on /3j through a tighter 
control of systematic uncertainties. 

As shown in Figure [26l these results are not strongly 
dependent on ^y. This is because the afterglow rate 
Faft is dominated by particles bouncing from the walls 
at grazing incidence, for which photon reflection itself 
contributes the large polarization-dependent phase shifts 
shown in Fig. [T] Averaging over polarizations weakens 
the dependence of Faft on <^y. 

The model-independent plots in this section and in 
Ref. [32[ assume that TOoS is the same in all runs. How- 
ever, in some models TOoff can depend significantly on the 
magnetic field. This is particularly true at the largest val- 
ues of Bq, where the electromagnetic energy density is of 
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FIG. 26: Model-independent scalar chameleon constraints for 
several reflection phases ^y. 
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27: Model-independent scalar chameleon constraints for 
of the magnetic fields used in CHASE, assuming = 0. 



1e+17 
1e+16 
1e+15 

_c 

Q. 1e+14 

o 

1e+13 



o 
o 
a. 



1e+12 



Collider constraints / / 



I 



1e+11 \r L^J. . 

Casimir force constraints 

1e+10 




ri,) 



x=io 

^=10"' 



1e+10 1e+12 1e+14 1e+16 
matter coupling 



1e+18 



1e+20 



FIG. 28: CHASE constraints on chameleons. Chameleon 
fragmentation becomes important for A > 0.01 and low /J^. 
Note that the entire parameter space shown has already been 
excluded by the Casimir force constraints of [4^ . 



low A « 0.01, fragmentation has only a negligible effect 
on afterglow constraints. Meanwhile, at low the con- 
straining power of CHASE is limited by the chameleon 
containment requirement; for low /3m and A, chameleons 
can escape though the walls of the chamber. 

As a result, CHASE constraints on 0^ chameleons 
are considerably weaker than fifth force constraints from 
measurements of the Casimir force (4^ . The entire pa- 
rameter space region shown in Fig.[25]is already excluded 
by Casimir force constraints. 



F. CHASE constraints on dark energy 



the same order of magnitude as the gas density inside the 
vacuum chamber. 

Thus, for completeness, we show in Figure [57] the con- 
straints resulting from each of the seven different mag- 
netic fields individually. A joint analysis of data and 
calibration runs is necessary to distinguish between the 
"orange glow" systematic and an actual chameleon after- 
glow, as described in Sec. lIXBl In order to avoid double- 
counting the calibration data, we analyze a different one 
of the seven calibration runs along with each magnetic 
field value in Fig. E71 



E. CHASE constraints on (p* theory 



Chameleon fragmentation in 0^ theory, V{(f)) ~ ji'P'^i 
was studied in Sec. IIII PI Figure ^ shows CHASE con- 
straints on (j)'^ theory. The effects of fragmentation are 
evident at large A and low /3-y, where fragmentation com- 
petes with oscillation as a means of reducing the de- 
tectable chameleon population in the experiment. Be- 



Next, we apply CHASE data to chameleon dark energy 
models ([9]). Since the constraints of the previous section 
apply also to n = 4 chameleon dark energy, we focus 
here on inverse power laws, n < 0. We also look at the 
exponential potential (j40p . which we write in a form 



Mi + Mi exp 



Kcf) 



(77) 



suitable to a dark energy model. As discussed in 
Sec. IIVBI these potentials are the easiest to study be- 
cause their phase changes s^tc independent of incident 
angle; more complicated potentials would require a nu- 
merical computation of as a function of 6. 

CHASE constraints on inverse power law chameleon 
dark energy are shown in Figure 1291 Twelve orders of 
magnitude in are excluded for a range of /3j in each 
of the models shown in Fig. [29l ( Top), with the excluded 
region limited at low by the containment require- 
ment and at high by destructive interference in os- 
cillation due to a large effective mass inside the chamber. 
Fig. [29] {Bottom) shows that changing k has the effect of 
translating the excluded region in the direction; the 
same behavior is seen for other n. Figure [SD] constrains 
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FIG. 29: Constraints on chameleon dark energy models ([9|. 
(Top) assumes k = 1 and varies n; (Bottom) fixes n — —1 and 
varies k. 
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FIG. 30: Constraints on exponential dark energy models (|77|l . 
Models with k > 10^ have particles too massive to be probed 
by CHASE. 



exponential dark energy (|77p . In these models nics = 
y/ ^(/SmPmat + PjPc,em) / (MaMpi) grows sorapidly with 
K that K > 10^ is inaccessible to CHASE 01. 
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FIG. 31: Combined constraints onn = — = 1 chameleon 
dark energy (|9|. Solid regions show current constraints from 
analyses of experimental data, while lines represent forecasts 
and preliminary analyses. 



G. Combined constraints on n = 1 dark energy 

Chameleon dark energy is described by a compli- 
cated parameter space, with several experimental con- 
straints which depend on the parameters in different 
ways. While the previous discussion attempted to com- 
pare CHASE to existing constraints, such constraints do 
not exist (or have not been published) for all of the po- 
tentials considered here. Furthermore, new experiments 
have constrained chameleon models since the publication 
of CHASE results in [s^ , and forecasts of planned exper- 
iments have been made. 

It is instructive to include as many chameleon con- 
straints as possible on a single plot. For that we choose 
a specific potential, chameleon dark energy (|9]) with 
n = — 1 and k = 1. Figure [21] shows current constraints 
and forecasts on this model. Collider constraints have 
been studied extensively [1^ H^, Hi], and are weakly 
dependent on /Sm and V{4>), but only exclude extremely 
strong photon couplings. Afterglow constraints are due 
to CHASE. At large matter couplings, the chameleon- 
mediated fifth force would affect electronic energy levels 
in atoms, which are inconsistent with the data [ij. More 
recently, fifth force contributions to the quantized energy 
levels of neutrons in a classical gravitational field have 
been used to constrain chameleons [45|, H^. Since the 
neutron wavefunction in such experiments is de-localized 
to a cloud of several microns, chameleon screening is in- 
effective for suppressing fifth forces, resulting in powerful 
constraints at large /3„i. 

The upcoming neutron experiments qBounce and 
GRANIT are expected to improve constraints on the 
matter coupling by several orders of magnitude (46l - l48j . 
"Helioscopes" designed to convert solar axions into pho- 
tons can also be used to constrain chameleons |49l - [5]| . 
Refs. (52l . [ssj study chameleon-photon oscillation using 
more distant astrophysical sources; however, since they 
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describe their analysis as "preliminary" and systematic 
effects due to astrophysical uncertainties remain to be 
analyzed, we include this constraint as a forecast. Not 
shown on the plot are Casimir force tests [s^l and the 
original GammeV Experiment |3l| . which do not con- 
strain the potential chosen here. The Eot-Wash torsion 
pendulum experiment [55| likely excludes a range of mod- 
els around = 1 , but constraints for this specific model 
have yet to be computed. 



X. CONCLUSION 

We conducted a thorough investigation of the physics 
of chameleon particles in afterglow experiments, focusing 
on questions which arose during the design and analysis 
of CHASE. Afterglow experiments rely upon two assump- 
tions about chameleon particles: reflection from dense 
matter, which is necessary for trapping chameleon parti- 
cles; and oscillation, which allows chameleon production 
as well as photon regeneration. By studying the interac- 
tion of a chameleon particle with atoms, in isolation as 
well as in a lattice, we showed that the matter density 
can be treated as homogeneous inside the chamber as well 
as in its walls. Chameleon-atom scattering is negligible 
in the vacuum; Fig. [5] showed that the chameleon parti- 
cle ignores the electron cloud, so that the cross section 
is dominated by hard-sphere scattering from the nucleus 
alone. Fig. [5] showed that a chameleon particle "sees" the 
chamber wall as a homogeneous solid until its mass be- 
comes much larger than the ^ 1 cV needed for reflection, 
so trapping is unaffected by the fact that real matter is 
a lattice of atoms. 

On the other hand, the smoothly varying magnetic 
field found in an afterglow experiment could poten- 
tially suppress chameleon production significantly. If the 
length scale on which the magnetic field drops from its 
maximum value to zero, at the edge of the magnetic 
field region, is much larger than the chameleon oscilla- 
tion length, then the transition will be adiabatic, and an 
incoming photon will emerge from the magnetic field in 
a pure photon state. This adiabatic suppression can be 
mitigated substantially by placing glass windows inside 
the magnetic field region, as is evident from a compari- 
son of Figs.lT^andfHl ( Top). Such interior windows were 
included in CHASE and will be essential to any future 
afterglow experiment which seeks to improve constraints 
at greater chameleon masses. 

After studying these general effects, we proceeded to 
calculate the signal expected in an afterglow experiment. 
The analytic approximation of js^ was improved upon in 



Sec.|Vl This approximation, in the limit where it becomes 
exact, was used to verify that the Monte Carlo simulation 
of Sec. IVIIll was correct. The Monte Carlo simulation, in 
turn, was used to show that the predicted signal will not 
be decreased by diffuse reflection in the chamber walls. 
The dependence of the signal on the wall reflectivity and 
the chamber geometry was shown in Fig. 1231 Afterglow 
and decay rates per chameleon particle for CHASE were 
shown in Fig. [21] 

Finally, in Sec. IIXI we explained in greater detail the 



analysis underlying the constraints of [32[. The model- 



independent constraints of that reference were shown for 
different chameleon phase shifts and for each of the mag- 
netic flcld runs used in the experiment. CHASE data 
were then used to constrain a wide range of models, in- 
cluding dark energy candidates. Since additional con- 
straints on chameleons have been released since the publi- 
cation of [12, and since further experiments are planned, 
we compared several different constraints and forecasts 
in Fig. [nH For the model shown, CHASE excluded 
five orders of magnitude in the photon coupling over a 
range of ten orders of magnitude in the matter coupling. 
CHASEhas made a substantial contribution to the study 
of photon-coupled chameleon field theories. 



Acknowledgments 

We are grateful to A. Baumbaugh, R. Cowsik, P. O. 
Mazur. R. Tomlin, A. Weltman, and W. Wester for 
many informative discussions. JS thanks the Brinson 
Foundation for its generous support. This work is sup- 
ported by the U.S. Department of Energy, Basic En- 
ergy Sciences, Office of Science, under contract No. DE- 
AC02-06CH11357. This work was supported by Fermi 
National Accelerator Laboratory, operated by the U.S. 
Department of Energy under contract No. DE-AC02- 
07CH11359. 

The submitted manuscript has been created by 
UChicago Argonne, LLC, Operator of Argonne National 
Laboratory (Argonne). Argonne, a U.S. Department of 
Energy Office of Science laboratory, is operated under 
Contract No. DE-AC02-06CH11357. The U.S. Govern- 
ment retains for itself, and others acting on its behalf, 
a paid-up nonexclusive, irrevocable worldwide license in 
said article to reproduce, prepare derivative works, dis- 
tribute copies to the public, and perform publicly and 
display publicly, by or on behalf of the Government. 



[1] E. Komatsu et al. Astrophys.J.SuppL, 192:18, 2011. ullah, et al. Astrophys. J., 746:85, 2012. 

[2] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. R. [4] Ariel G. Sanchez, C.G. Scoccola, A.J. Ross, W. Percival, 

Nolta, et al. Astrophys. J. SuppL, 192:16, 2011. M. Manera, et al. 2012. e-Print arXiv: 1203.6616. 

[3] N. Suzuki, D. Rubin, C. Lidman, G. Aldering, R. Aman- [5] L. F. Abbott. Phys. Lett. B, 150:427, 1985. 



29 



[6] J. D. Brown and C. Teitelboim. Phys.Lett, B195:177- 

182, 1987. [32 
R. Bousso and J. Polchinski. JHEP, 6:6, 2000. 
P. J. Steinhardt and N. Turok. Science, 312:1180-1182, [33 
2006. 

C. de Rham et al. Phys. Rev. Lett, 100:251603, 2008. [34 

C. de Rham, S. Hofmann, J. Khoury, and A. Tolley. 
JCAP, 0802:011, 2008. [35 

G. Dvali, S. Hofmann, and J. Khoury. Phys. Rev., [36 
D76:084006, 2007. 

N. Agarwal, R. Bean, J. Khoury, and M. Trodden. [37 
Phys. Rev., D81:084020, 2010. 

E.G. Adelberger et al. Prog.Part.Nucl.Phys., 62:102-134, 
2009. [38 
J. A. Frieman, C. T. Hill, A. Stebbins, and I. Waga. Phys. 
Rev. Lett, 75:2077-2080, 1995. [39 
S. M. Carroll. Phys. Rev. Lett, 81:3067-3070, 1998. 
J. Khoury and A. Wehman. Phys. Rev. Lett, 93, 2004. [40 
171104. 

J. Khoury and A. Weltman. Phys. Rev. D, 69, 2004. [41 
044026. 

Ph. Brax, C. van de Bruck, A.-C. Davis, J. Khoury, and [42 
A. Wehman. Phys. Rev. D, 70, 2004. 123518. 
S. S. Gubser and J. Khoury. Phys. Rev. D, 70, 2004. 
104001. [43 
A. Upadhye, S. S. Gubser, and J. Khoury. Phys. Rev. D, [44 
74, 2006. 104024. [45^ 

D. F. Mota and D. J. Shaw. Phys. Rev. Lett, 97, 2006. [46 
151102. 

D. F. Mota and D. J. Shaw. Phys. Rev. D., 75, 2007. [47 
063501. [ 
K. Hinterbichler and J. Khoury. Phys. Rev. Lett., 
104:231301, 2010. [49 
K. A. Olive and M. Pospelov. Phys. Rev., D77:043524, [50 
2008. 

Massimo Pietroni. Phys. Rev., D72:043535, 2005. [51 
C. Deffayet et al. Phys. Rev., D65:044026, 2002. 
A. Nicolis, R. Rattazzi, and E. Trincherini. Phys. Rev., [52 
D79:064036, 2009. 

P. Brax, C. Burrage, A.-C. Davis, D. Seery, and [53 
A. Wehman. Phys. Rev. D, 81:103524, 2010. e-Print 
arXiv:0911.1267. [54 

H. Gies, D. F. Mota, and D. J. Shaw. Phys. Rev. D, 
77:025016, 2008. 

M. Ahlers et al. Phys. Rev. D., 77:015018, 2008. [55 
A. S. Chou, W. C. Wester, A. Baumbaugh, H. R. 
Gustafson, Y. Irizarry-Valle, P. O. Mazur, J. H. Stef- 
fen, R. Tomlin, A. Upadhye, A. Weltman, X. Yang, and 



J. Yoo. Phys. Rev. Lett, 102, 2009. 030402. 

J. H. Steffen et al. Phys. Rev. Lett, 105:261803, 2010. 

ePrint: arXiv: 1010.0988. 

J. H. Steffen and A. Upadhye. Mod. Phys. Lett. A, 
24:2053-2068, 2009. 

A. Upadhye, J. H. Steffen, and A. Weltman. 

Phys. Rev. D, 81:015013, 2010. 

J. H. Steffen et al., 2012. (in prep.). 

G. Raffeh and L. Stodolsky. Phys. Rev. D, 37:1237-1249, 

1988. 

P. Brax, C. van de Bruck, A. C. Davis, D. F. Mota, 
and D. J. Shaw. Phys. Rev. D, 76:085010, 2007. e-Print 
arXiv:0707.2801. 

W. A. Rolke, A. M. Lopez, and J. Conrad. Nucl. In- 
strum. Meth. A, 551:493-503, 2005. 

K. Hinterbichler, J. Khoury, and H. Nastase. JHEP, 
1103(61), 2011. 

P. Brax, C. Burrage, A.-C. Davis, D. Seery, and A. Welt- 
man. JHEP, 0909:128, 2009. e-print arXiv:0904.3002. 
M. Kleban and R. Rabadan. 2005. ePrint arXiv:hep- 
ph/0510183. 

P. Brax, C. van de Bruck, A. C. Davis, D. F. Mota, 
and D. J. Shaw. Phys. Rev. D, 76:124034, 2007. e-Print 
arXiv:0709.2075. 

R. Balest et al. Phys. Rev. D, 51:2053, 1995. 

P. Brax and C. Burrage. Phys. Rev. D, 83:035020, 2011. 

V. V. Nesvizhevsky et al. Nature, 415:297, 2002. 

P. Brax and G. Pignol. Phys. Rev. Lett, 107:111301, 

2011. 

M. Kreuz et al. 2009. ePrint arXiv:0902.0156. 

S. Bassler, A. Gagarski, L. Grigorieva, M. Kreuz, 

F. Naraghi, et al. 2012. e-Print arXiv: 1202.2784. 

P. Brax and K. Zioutas. Phys. Rev. D, 82:043007, 2010. 

P. Brax, A. Lindner, and K. Zioutas. 2011. ePrint 

arXiv:1110.2583. 

O. K. Baker, A. Lindner, A. Upadhye, and K. Zioutas. 
2012. ePrint arXiv: 120 1.0079. 

C. Burrage, A.-C. Davis, and D. J. Shaw. Phys. Rev. D, 
79:044028, 2009. 

C. Burrage, A.-C. Davis, and D. J. Shaw. 
Phys. Rev. Lett, 102:201101, 2009. 

Ph. Brax, C. van de Bruck, and A. C. Davis. 
Phys. Rev. Lett, 99, 2007. 121103 [arXiv:hep- 
ph/0703243]. 

D. Kapner et al. Phys. Rev. Lett, 98:021101, 2007. e- 
Print arXiv:hep-ph/0611184. 



